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The limitation on the Geiger-Miiller counter insensitive time which is imposed by the presence 


of the positive ion space charge has been reduced by collecting the positive ions on the counter 
center wire instead of at the outer cylinder. The electronic circuit which carries out this process 
is described in detail. A reduction by another order of magnitude in the insensitive time or an 
insensitive time of about 2X 10~ sec. is obtained with no indication that the limit of the new 
method has been reached. The general properties of an inverse counter are discussed. The 
mechanism of this mode of counter operation is explained with considerable evidence indicating 
the presence of three distinct ion collection regions for the negative center wire. A stable circuit 


useful as a research tool has been devised for use with a large variety of counters. 


I. INTRODUCTION 


AREFUL studies have been made'? of the 
positive ion sheath which show that there is 
a natural insensitive or dead time associated with 
each type of G-M counter. This insensitive time 
exists because the positive ion space charge re- 
duces the electric field in the center-wire region 
below the critical value for supporting a new 
avalanche until the space charge has moved 
outward about two-thirds of the distance toward 
the cathode. If a high mobility ion sheath could 
be formed, this transit time would be reduced 
considerably. Such a sheath is encountered in a 
permanent gas mixture of two light gases, but to 
date no method has been successful in preventing 
secondary electrons and photons from forming in 
such a counter. If an attempt is made to use 
polyatomic gases which have a very small mass 
and size, then the dissociation process from 
* Now at the University of Chicago, Chicago, Illinois. 
'C.G. Montgomery and D. D. Montgomery, Phys. Rev. 


57, 1030 (1940). 
*H. G. Stever, Phys. Rev. 61, 38 (1942). 


counting results in continual changes with time 
of the counter characteristics, and the counter has 
a very short lifetime. 

This insensitive time of G-M counters places a 
severe limitation on their use in obtaining accu- 
rate data at high counting rates. Since, in general, 
the events which are recorded by such counters 
are purely random, at high counting rates there 
will be a large number of intervals between two 
successive events which are less than the insensi- 
tive time interval of any counter. This loss of 
counts at high rates is expressed in an equation 
derived by Ruark and Brammer* in which the 
insensitive time is o; the true counting rate is NV, 
and the observed rate is N,; the equation is 


In order to obtain some idea of these losses for 
different values of o, the family of curves for Na, 
(N./N.)—1 with o as a parameter appears in 
Fig. 1. Some values for ¢ the order of 10-* second 


+A. E. Ruark and F. E. Brammer, Phys. Rev. 52, 322 
(1937). 
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have been chosen which agree with experimental 
values obtained for so-called ‘‘fast’’ counters. 
The very best conditions that can be obtained in 
recording data from a counter with insensitive 
time ¢ are attained when the amplifier responds 
to the very tiny pulses which may occur im- 
mediately after the elapse of the time interval ¢ 
following each ionizing event, unless the counter 
‘is operated just below the Geiger threshold. In 
the latter case, the sheath does not fully develop, 
and the pulses are not independent of the kind of 
ionizing event. 

A solution to the problem is usually not found 
by changing the geometry of the sources of the 
ionizing events during the experiment, especially 
if the high counting rate is changing rapidly with 
time. If the source is placed at such a distance 
that it produces a counting rate low enough to 
avoid serious losses, then when the source decays 
rapidly, the counting rate will become so reduced 
that exceedingly long times are required to main- 
tain the accuracy of the individual measurements. 


Il. A NEW COUNTER MECHANISM 


Another mode of operation of a G-M counter 
is evidently necessary if progress is to be made in 
reducing o. It has been effectively shown by 
Ramsey‘ that the total collection time for elec- 
trons on the center wire of an organic vapor-noble 
gas counter is less than 10~* second, and most of 
the electrons are collected in about this time in a 
“‘slow’’ counter. Therefore, after the center wire 
begins to drop in potential, use is not made of the 
subsequent cycle so far as recording the event is 
concerned. Consider an alternative procedure. 
Instead of allowing the positive ion sheath which 
forms around the center wire to move outward, 
suppose these ions were almost immediately col- 
lected on the center wire. This may be accom- 
plished by reversing the potentials on the center 
wire and cylinder long enough to insure collection. 
Since the ions are formed only in the immediate 
vicinity of the wire and will therefore be in a high 
field which accelerates them toward the center 
wire, this collection time will be very short, 
probably less than 10-° second. As soon as all 
ions are collected, the counter should be im- 
mediately returned to its normal operating po- 


*W. Ramsey, Phys. Rev. 57, 1022-1029 (1940). 
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tential to be ready for the next event. In this case 
the outer cylinder serves only as a field-forming 
conductor. 

It is important to consider the properties and 
action of a counter when it is inverted. As in the 
case of the positive center wire, the negative 
center-wire counter has a definite potential at 
which it breaks into corona. For the impure gases 
which exist in G-M counters, the breakdown 
potential at very low pressures when the center 
wire is positive is greater than when the center 
wire is negative. The pressure-potential curves, 
however, cross at higher pressures, and for pres- 
sures used in most counters the breakdown po- 


Actuol Per Sec 


Fic. 1. This is a graph to determine the true counting 
rate N;, from the observed counting rate N, for a given 
insensitive time 


tential of the negative center-wire system will, in 
general, be greater. 

Below the breakdown potential in a negative 
center-wire system there is a small region in which 
the system behaves like a G-M counter in some 
respects. An ionizing event in the volume results 


- generally in a single pulse, and the pulse is inde- 


pendent of the number of primary ions produced. 
This may be called the inverse Geiger region. 
Counters operated in this fashion have been tried 
principally by Ramsey and Cowie.* However, it 


5M. J. Druyvesteyn and F. M. Penning, Rev. Mod. 
Phys. 12, 122 (1940). These authors present a complete 
brief discussion on breakdown and corona with a concentric 
electrode system. 

*D. B. Cowie, Phys. Rev. 48, 883 (1935). 
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is evident that most ionizing events take place in 
the outer volume of the counter, and therefore 
the time before a positive ion reaches the center- 
wire region is a variable; the positive ion hasa 
very low mobility compared to electronic mo- 
bility. As a result, negative wire counters cannot 
be used in a coincidence circuit. The background 
of such a counter is high, and it appears to have 
no appreciable plateau. If the potential is raised 
slightly, but not enough to cause breakdown, the 
counter may break into periodic groups of high 
frequency oscillations. The envelope of the 
groups is well defined. An explanation of this will 
be discussed below. 

By making the center wire negative, the pro- 
duction of photoelectrons becomes negligible. 
Photoelectrons are produced at the cathode of the 
electrode system which in this case is the center 
wire. Therefore, this photo-effect is reduced by a 
factor approximately equal to the ratio of the 
wire radius to the cylinder radius. It should be 
noted that tungsten, which is widely used for 
center wires, has a work function of about 4.5 ev. 

The behavior of positive ions near the negative 
center wire for potentials below the inverse 
Geiger region is probably quite complex and does 
not lend itself to an exact solution. It is desirable 
to know under what conditions the positive ions 
may acquire enough energy to produce secondary 
electrons by field emission at the wire. If the ions 
are polyatgmic, they may be polarized near the 
wire and change the effective mean-free-path in 
that region. At the conclusion of this discussion 
something will be said concerning these processes 
in the light of experiments to be presented. 

The avalanche process at the negative wire has 
many characteristics which differ from the well- 
known positive center-wire avalanche. If a posi- 
tive ion approaches the center wire and causes 
the ejection of an electron by field emission, this 
electron will travel toward the outer cylinder. 
The field is decreasing as 1/r. By the time electron 
avalanche production would be able to increase 
to very large values, as in an increasing field, the 


_ avalanche head has already entered a low in- 


tensity field region. Townsend’s coefficient a, 
which is a measure of the number of ion pairs 
formed by an electron in one centimeter of its 
path, depends upon the field intensity and the 


INSENSITIVE TIME IN G-M COUNTERS 


41 


pressure. Since a@ is decreasing almost expo- 
nentially as the avalanche head moves into a 
region of low E/p, the production will break off 
fairly sharply. The field is further distorted and 
reduced by the positive space charge left behind . 
the avalanche. The current flow of electrons, 
therefore, increases rapidly to a value which is 
constant for most of the distance out to the 
cylinder. If many electrons initiate a large num- 
ber of avalanches, the positive space charge may 
be sufficient momentarily to choke off further 
avalanche production until a large portion of the 
positive ions have been collected; the process 
then repeats. This would produce an oscillatory 
discharge in the counter of the nature previously 
described. 

The negative ions play an important part if a 
great number exist when the counter potentials 
are reversed, for then the negative ions must 
travel out toward the cylinder. However, if only 
a few such ions are formed, they need cause no 
difficulty because they will not distort the field 
distribution ; they will not lose their electrons in 
a region where a new avalanche could possibly 
start. The field in front of the outer cylinder is 
very low. 

These general considerations indicate that the 
new mode of counter operation proposed can be 
successful if the proper inverted region is se- 
lected. Later, this will be shown to be possible. 


Ill. THE EXPERIMENTAL CIRCUIT 


The requirements for the design of an electronic 
circuit to produce the random reversals at a very 
high average rate with precision are: 


1. The input to any amplifying device which 
operates on the first three or four volts’ drop of 
the steep pulse front must be insensitive to the 
process of reversing the field. 

2. The circuit must have extremely low RC 
times throughout; it must respond to frequency 
components around two megacycles and still 
possess fairly high gain. 

3. The reversal process must be carried out so 
that not more than a few microseconds elapse 
before the field is inverted and ions are being 
collected. 

4. No residual voltage ‘‘tails’’ more than three 
or four volts can be tolerated after collection time 


} 
ase 
ing 
nd 
he 
ve 
at 
wn 
ter 
ter 
es, 
ig 
n 
n 
h 
e 


42 J. A. SIMPSON, JR. 


is over or, at high counting rates, the effective 
operating potential will be changed. 

5. Accurate control of the time and extent of 
the reversal of field must be possible over a wide 
range to accommodate all types of counters. 

6. All supply potentials must be constant 
under extreme fluctuations of counting rates. 

7. The circuit should be constructed from 
standard equipment available to all laboratories 
if it is to be of use as a research tool. 


A circuit fulfilling these requirements has been 
developed and operated successfully ; before de- 
scribing it in detail, it would be best to outline the 
general method with the aid of a simplified block 
diagram, Fig. 2, representing the circuit elements. 
Normally no current flows in R,; therefore the 
counter is operating at a potential V2. When the 
electrons are collected from an ionizing event, a 
positive pulse begins to form across R. This is 
detected, and a sharp pulse passes through to the 
first trigger pair which produces a rectangular 

‘wave of time length ¢;; this is fed back to the 
input circuit to act as a shutter preventing the 
first portion of the circuit from functioning for a 
time f;. The second trigger pair is actuated by the 
incoming pulse also, and it produces a rectangular 
wave of width ¢:. This rectangular pulse is con- 
verted to a rectangular pulse approximately 
—(Vi+V:2) volts high appearing at a. If >t, 
then all disturbances on the center wire are over 
before the shutter opens again. 

The schematic diagram of the circuit actually 
used appears in Fig. 3. All values of the circuit 
constants are listed. A beam power tetrode of 
type 807 (or type 829 with its elements in 

parallel) is used for the high voltage pulse- 
forming tube with the grid bias of about —150 
volts adjusted to place the grid potential beyond 
cut-off ; thus the average plate current required is 
very low. With the center wire connected directly 
to the plate of the tube, the full counter operating 
potential V2 is applied between ground and plate 
supply potential. The cathode and heater are 
— V; volts below ground potential with Vi+ V2 
volts applied across the type 807 tube. The series 
of pulses to be described appears in Fig. 4A. The 
positive pulse front a which is several volts high 
in one or two microseconds appears across the 
resistor R since the impedance in the plate circuit 


of the 807 is negligible in comparison with R. 
Pentode A is biased just beyond cut-off, so it will 
produce a sharp negative pulse b across its plate 
impedance when the small positive pulse (1-3 
volts) is applied directly to the grid of A. 
Pentode B amplifies and inverts this pulse c by 
operating well up on its normal characteristic 
curve. From here, the positive pulse is applied to 
the tube C which serves two purposes. It produces 
a negative pulse d for the trigger pair D and E; 
it also is used as the shutter for the amplifier by 
receiving a rectangular wave e from the trigger 
pair on the biased suppressor grid. The shutter 
time ¢, is controlled by Ris. The trigger pair 
operates in the conventional manner with tube D 
normally conducting and with E normally non- 
conducting. The steep front of the negative rec- 
tangular pulse is differentiated f by the circuit 
RisCy. This triggers the tubes F and G which are 
arranged exactly as D and E so as to produce a 
rectangular pulse g whose time width fz is de- 
termined by Ris. The pulse is now amplified by H 
and inverted h so it can drive the grid of the type 
807 so far positive that small variations in the 
height of A will not affect the characteristics of * 
the negative pulse i of approximate height 
volts. But the outer cylinder of the Geiger counter 
is near ground potential ; therefore the center wire 
drops to approximately V, volts below ground 
potential. Since the plate circuit has a much 
lower impedance than R, the differentiated rec- 
tangular wave will appear at the input of the 
amplifier tube A, but a time ¢;>¢: must elapse 
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Fic. 2. This block diagram illustrates the manner in 
which the circuit elements are arranged to produce positive 
ion collection on the Geiger-Miiller counter center wire in 
order to reduce the counter insensitive time. 
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Fic. 3. Schematic diagram of electronic circuit to collect positive ions on the center wire after each ionizing event. The 


circuit constants are as follows: 


C, =0.1 uf, Cro =0.1 yf, 

Cu = 100 pul, 
= 0.02 uf, 
Cis = 500 
Cus =0.001 uf, 
Cu =4 uf, 

Ly =1.2 mh, 
L. =0.67 mh 
L; =1,.2 mh, 


L, =80 mh, 
R, = 35000, 
R. = 5002, 

R; =25 K, 

R, =25 K, 

R; = 25002, 
R, =25 K, 

R; =0.1 meg., 
R, =25 K, 


before the amplifier passes any pulses through the 
shutter, and so this pulse is unrecorded. 

A vertical section through the diagram (Fig. 5) 
of the potential changes on the counter electrodes 
as a function of time gives the potential between 
the electrodes at any instant. The shaded area is 
drawn between the two electrode potential 
curves. The first pulse a is purposely drawn 
incorrectly to point out the importance of avoid- 
ing ‘‘tails’’ on the pulses. With pulse a the counter 
does not get back onto its operating potential 
until the tail of the reversing pulse has disap- 
peared. The tiny differentiated pulse appearing 
at the outer cylinder further tends to reduce the 
operating potential immediately after ion collec- 
tion. Pulse b is the composite of the two pulses 
actually observed on the oscilloscope when the 
circuit is functioning correctly. 


IV. CIRCUIT TESTS 


The preliminary testing of circuit elements is 
of importance. The tubes A, B, and C were over- 


Ru=4 K, 
Res = 25000, 
Rx» = 5002, 
R =75,000. 


Ry =2.5 K, 
Ryw= 25 kK, 

Ry =75 K w.ew. 
K, 
Ri3=2 meg., 
Ryu= 10 K, 
Ris=4 kK, 

50 K, 
Ru=6 K, 


Ris=1 meg., 
Riy=2 K, 
50 K, 
Run =6 kK, 
meg., 
Ry3= 5002, 
Ry =0.1 meg., 
K, 
Ro = 50 


A, B, C, D, E, F 
G, Hare 6AC7. 


compensated with inductance to avoid severe 
losses of the higher frequency components of the 
pulses by tube and circuit capacities. Such ampli- 
fiers must be untuned because of the aperiodic 
nature of the pulses. The upper frequency limit 
of the amplifier is many times that of the pulses 
themselves. The very low frequencies were of no 
particular concern in this work. Decoupling in 
the plate circuits is provided wherever it might 
possibly avoid difficulties. Tests were made with 
microsecond pulses over a wide frequency range. 
The power supply for tubes A through H was a 
single unit run from a stabilized power line. 
The shutter circuit composed of tubes C, D, 
and E must be carefully adjusted. The biases on 
the control grid and suppressor grid are so ad- 
justed that the tube is normally just beyond cut- 
off. A pulse appearing on the control grid is 
amplified, but a negative rectangular wave ap- 
plied to the suppressor grid will keep the tube 
non-conducting. This circuit element is checked 
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by applying sharp microsecond pulses at about 
100-200 kc to the control grid while varying the 
time ¢;. By so doing the frequency can be divided 
by the shutter by any integral factor up to the 
limit of ¢, determined by tubes D and E. The 
input to D must be carefully adjusted. Several 
methods of coupling the pulses from F and G to 
the type 807 have been tried with success, but the 
one illustrated was used in obtaining the results 
described in this paper. 

Since the grid is Vi;+150 volts below ground, 
it is essential that C,, be a high voltage condenser 
and that the separate filament transformer be 
insulated if widely varying potentials of V, are 
desired ; otherwise, for normal operation no such 
precautions are necessary. The screen supply 
voltage may conveniently range from 300-600 v. 
On counters with low operating potentials a 6C6 


tube or equivalent type has worked with some 


success. 

It has already been pointed out that no tail 
(>3 or 4 v) must be present on the negative 1000 
to 1500 v pulses produced by the type 807 tube. 
If only a plate resistor with the usual small 
compensating inductance were used, long tails 
starting about 150-200 v in height appear even 
with small values of R,. This problem was 
effectively solved by means of an LRC circuit. 
Used with a variable inductance L, (approxi- 
mately 15 millihenries), the resistance R is 


Fic. 4A (Above). The pulse shapes observed on an 
oscilloscope connected to different portions of the circuit. 
The pulses are not drawn to scale. The input pulse is the 
first positive front of a; the reversing pulse is 1. 

Fic. 4B (Below). The effect of varying the value of the 
inductance L, and its shunt resistance R, on the tail of the 
reversing pulse. The tail must be very small so that the 
counter can return to its operating potential rapidly. 


Fic. 5. A graphic illustration which indicates the 
potential difference between the center wire and outer 
cylinder of the counter at any instant of the counter 
reversal process. The shaded areas are drawn between the 
two electrode potential curves. 


critical and has a value of about 550 ohms. 
Figure 4B shows the results (not drawn to scale) 
of varying both the shunting resistance and the 
inductance. The tail in pulse VI is less than 4 v 
high above the axis and disappears within the 
time 

The separate power supply for the type 807 is 
stabilized electronically. Under a load of 100 to 
150 milliamperes the voltage should not fluctuate. 

By using periodic groups each of two pulses of 
width 10-6 second, spaced 10~* second apart, the 
circuit is checked for selectivity. The second 
pulse is eliminated, and the circuit only reverses 
on the first pulse of each group as it should. The 
whole circuit has been tested at frequencies over 
300 ke with a coupling condenser instead of a 
counter between tube 807 and the input of tube A 


_ which has the same value as the counter capacity. 


All observations were- made on a Dumont 
oscilloscope No. 241, a shielded probe, and high 
impedance attenuator. Care must be taken that 
the observation does not alter the circuit in any 
manner. Any oscilloscope whose amplifier charac- 
teristics are not flat in frequency response up to 
2.5 megacycles should not be used. 

The scaling and recording circuit is operated 
from the output at E. The first scale of 8 or 16 
functions up to about 3 megacycles, and the re- 
maining scalers to form a scale of 512 are con- 
ventional hard-tube circuits. 

To find the time delay in making the center 
wire V; volts negative with respect to ground, the 
sharp input pulses were superposed with the re- 
versing pulses on the oscilloscope screen. This 
delay is about 1.5 X 10~* second for a 1000-v pulse. 
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V. RESULTS 


The reduction of the insensitive time of a G-M 
counter by this circuit can be determined by any 
one of several methods for measuring resolving 
times. The most useful one for instantaneous ob- 
servation is the single-sweep trigger circuit as 
used by Stever.? One event triggers off an oscillo- 
scope to sweep across the screen ; at high counting 
rates there is a fair probability that the next 
event will occur near the limit of the insensitive 
time. This is observed by pulses making their ap- 
pearance on the screen after the insensitive time 
interval. Thus, by repetition of the pattern the 
interval between the trigger pulse and the next 
appearance of the pulses measures the insensitive 
time o. By using this circuit connected to the 
first or second trigger pair in the reversing circuit, 
the circuit characteristics will not be changed. 
All pulses will be of the same height in contrast to 
the usual observations.? These observations have 
been made with a calibrating time scale (3 per- 
cent error) incorporated in the sweep. 

Polyatomic-noble gas counters are most useful 
for the first experiments for two reasons. When- 
ever the reversing circuit is shut off, the counter 
will not be destroyed. Also this type of counter 
has the lowest range of insensitive times used at 
present, and therefore it is a suitable standard for 
comparison with results obtained by using the 
reversing Circuit. 

The factor of reduction of insensitive time is 
plotted in Fig. 6 as a function of the magnitude 
of the reversing potential for the first counters’ 
used in the circuit. ¢_or is the insensitive time of 
the counter without the circuit; o_o, is the 
insensitive time of the counter when the reversing 
circuit is used. Great care was taken that the 
smallest pulses would be recorded. For each 
counter V2 isaconstant. Preliminary experiments 
indicate that the new insensitive times lie around 
2X 10-5 second. This is nowhere near the limit to 
be expected in future trials. It should be noted 
that these values are obtained under very stable 
conditions ; the apparatus may be turned off and 
warmed up again as desired without difficulty. 
Background counting rates remain a constant in 
using the circuit. 

7 Counter No. 2 has 10 percent amyl acetate; 90 percent 


argon. Counter No. 5 has 10 percent ethyl alcohol; 90 
percent argon. - 


In obtaining the points for the curves in Fig. 6 
the collection time was adjusted for its optimum 
value for each value of V; with the shutter pulse 
set at 4;=1.510-5 second. The time ¢2 varied 
from 2.5X10-* to about 6X10-* for the ion 
collection. 

Additional tests were made by using two y-ray 
sources, one twice the intensity of the other. The 
intensities of these sources are first accurately 
determined. Then the lower intensity source 
(No. 1; 9.98 mg) is determined at, say, 10 counts 
per second ; the second source (No. 2; 19.98 mg) 
is placed in the same position, and the counting 
rate is again determined. Finally, with both to- 
gether the counting rate is obtained. After re- 
peating this procedure at several distances from 
the counter the information can be used to obtain 
curves like those in Fig. 1. 

By varying the potential and width of the re- 
versal pulse three distinct regions in the inverse 
Geiger region are found in polyatomic-noble gas 
counters. This has been anticipated to some extent 
in the introductory remarks of this discussion. 
At low voltages (up to about 300 v with counters 
operating at around 1100 v) the ion sheath is 
collected within about 7X10-* second without 
producing any detected secondary effect. The 
upper limit of this potential is defined within 
about 40 v, but its magnitude depends upon the 
counter properties, especially the mean free path 
at the center wire. The diameter of the center 


BH * * secon 


as 

Fic. 6. Experimental evidence of the decrease in counter 
insensitive time produced by reversing counter electrode 
potentials for positive ion collection. ¢o¢¢ is the insensitive 
time with the reversing circuit off; oo, is the insensitive 
time using the circuit. V: is the counter operating potential ; 
V, is approximately equal to the reversed potential. 
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wire will also determine the terminus of the first 
region. An experiment in which the diameter of 
the center wire is varied has been planned. 

The second region begins when tiny sharp 
pulses (less than 0.3 v) appear at the outer 
cylinder about 2X 10~‘ second after the ionizing 
event a. Figure 7 illustrates the changing charac- 
ter of this ‘‘tail” as a function of Vi, which ap- 
pears at the end of the differentiated negative 
pulse. In each successive drawing the potential V; 
has been increased. It is possible to hazard an 
interpretation of this phenomenon. In a and 6 the 


small pulses are due to positive ions arriving at 


the outer cylinder. The beginning of the second 
region is one in which the noble gas positive ions 
begin to have enough energy to produce second- 
ary electrons at the center wire; these electrons 
will produce new local avalanches whenever the 
field is intense enough to support them. The space 
charge is being collected before transformation 
by electron transfer can occur; therefore a con- 
siderable portion of the sheath is composed of the 
noble gas ions. An example of this process appears 
in the literature. Penning® has observed that the 
number of electrons produced per positive ion 
impact varies from about 0.03 to 0.05 for ions 
having energies less than 50 v. This fraction in- 
creases nearly linearly with energy up to 0.42 for 
1000-v ions. 

How are the characteristic patterns of a, 6, and 
c to be explained in view of the fact that the ion 
collection times can be varied somewhat without 
changing them? As the positive ion sheath is 
coming toward the center wire, the field between 
sheath and wire becomes more intense, but the 
region beyond the sheath drops to a relatively 
low value because of this space charge. Secondary 
electrons passing outward are, therefore, unable 
to start avalanches in this region. As the sheath 
is being collected, the field again rises to the point 
where local avalanches may be produced. The 
continued collection of the sheath brings new 
local avalanches closer and closer to the wire. If 
the counter were suddenly returned to normal 
conditions at this instant, small individual space 
charges would move out to the cylinder to be 
collected, thus producing small peaks in the ion 
collection pulse. But if the counter were not re- 


8F. M. Penning, K. Akad. Wet. Amst. Proc. 33, 841 
(1930). 
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versed at this instant, the new positive ions 
would arrive at the wire to continue the process 
until the counter potentials were changed. Again 
the ion collection pulse would be a series of 
peaks. Little or no spread of the discharge due to 
photons may be expected in the lower part of this 
region in which only a few of the positive gas ions 
have enough energy to produce secondary elec- 
trons. The resolution of pulses in a and b tends to 
confirm the assumption made in G-M counter 
theory that the ions move out without much 
spreading, i.e., the sheath can be treated as 
infinitely thin for purposes of calculation. 

As the potential V; is increased, secondary 
electrons are produced all along the wire while 
positive ions are being more rapidly removed 
from the original sheath. A large space charge 
from the electron avalanches then forms not only 


Fic. 7. Sketches of the pulses observed by using an 
oscilloscope connected to the outer cylinder. The potential 
V; is great enough to cause positive ions to produce 
secondary electrons at the center wire. V; is increasing from 
left to right. 


around the wire, but also along its length; the 
role of photons in this process is questionable. 
The large space charge immediately chokes off 
further avalanche production, and the large burst 
of electrons which formed the avalanche shortly 
arrives at the outer cylinder. The collection of 
these electrons can be observed by the negative 
pulse appearing at a in d, e, and f. The magnitude 
of this pulse is a function of the potential V; and 
not of the counter overvoltage; therefore it can- 
not be interpreted as a collection of the original 
positiye ions. This pulse is terminated within 
4-5X10-* second from the time of reversal. It 
demonstrates that the original positive ion sheath 
is at least mostly collected within that time 
interval. When the time fz is extended, a repeti- 
tion of the process takes place and a second 
negative peak is to be observed: When the 
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counter potentials are returned to normal oper- 
ating potential, this extensive positive ion sheath 
moves outward producing the positive hump 
shown on the tail of the differentiated pulse. A 
counter operated in such a region obviously has 


its insensitive time increased rather than de- ° 


creased. 

The negative pulse a and the positive hump 
continue to increase with the potential of V; up 
to a value | Vi|=2| V2| where the third region 
begins. This is the inverse Geiger region described 
in the introduction. Multiple and spurious pulses 
make their appearance. Above this potential lies 
the corona discharge region which is well known 
and is useless in these observations. 

By triggering the reversing circuit artificially 
in any of these regions at a time when no ionizing 
event is present, it can be shown that the differ- 
entiated pulse lacks all the characteristic pulses 
described in the preceding paragraphs; the phe- 
nomena are due to counter action and not to the 
characteristics of the reversing pulse. In Fig. 6 
the point where the curve crosses the unity axis 
indicates the vicinity of the boundary between 
the ion collection region and the region in which 
ions are beginning to produce secondary electrons. 

Another way in which to collect the positive 
ions at high counting rates may be pointed out. 
Suppose the circuit reverses periodically, being 
inverted for a time ¢’ and returning to normal for 
a much longer time ¢’’. Further, ¢’’ will be made 
smaller than the normal insensitive time tz. A 
complete cycle of operation, therefore, requires a 
time ¢’+?¢’<tg. Conservative values for these 
times can be chosen, such as 


t'=10-* second, 
second. 


To a first approximation the probability of two 
events in a random distribution occurring in a 
time interval less than + with a counting rate 
N/sec.=1/t is* p=7/t. Assuming a counting 
rate of 500 per sec. the probability of two events 
occurring in a time interval sec. is 510%, 
in a time interval 5X10-5 is 2.5X10-, and for 


*]. Strong, Procedures in Experimental Physics (Prentice- 
Hall, Inc., New York, 1938). 


10-* sec. (the dead time of a “fast’’ normal 
counter) is 5X10-* sec. The advantage of such a 
procedure over using the counter alone is further 
emphasized since it is definitely known that } of 
all the counts are purposely not being recorded. 

Although the circuit in Fig. 3 was designed to 
study the properties and possibilities of reducing 
the insensitive time, it may be used as a research 
tool ; there are several critical adjustments which 
have been described, but once adjusted, the 
circuit is stable. For use with a particular counter 
type a much simpler circuit may be devised. The 
gain requirement may be relaxed. A simple 
design would include two stages of amplification 
incorporating the shutter along with two LC cir- 
cuits to produce rectangular waves of the desired 
widths and a pulse amplifier followed by a 
blocking oscillator using a pulse transformer. 
The large tail produced by such an oscillator may 
be corrected by the LRC circuit suggested earlier. 
The use of a high voltage trigger pair to obtain 
the reversing pulse has the disadvantages of re- 
quiring a large average plate current and pro- 
ducing pulses of an undesirable shape. The new 
techniques which will be available in the future 
for pulse production may result in a drastic 
simplification of the reversal circuit. Research is 
still to be carried out on the permanent gas 
counter ; the advantage of using such a counter is 
that it should maintain its characteristics until 
the wire becomes badly pitted. Since the outer 
cylinder now only assumes the character of a 
field-forming electrode, serious consideration 
should be given to the fact that almost any sur- 
face layer which does not accumulate charges 
may be used to make the production of photo- 
electrons impossible. 
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The energy spectrum of electrons at the altitude of 14,618 feet originating from the decay of 


mesotrens is calculated by means of the cascade theory of showers. It is assumed that the 
mesotron decays into a neutrino and an electron, which are emitted with equal probability in 
all directions in the reference frame of the mesotron. The diffusion equations for the showers 
are derived with the boundary conditions that there are no incoming electrons or photons. 
By approximations similar to those used in the shower theory, it is possible to calculate the 
electron spectrum in the lower atmosphere, provided the mesotron spectrum is known. If we 
assume a mesotron spectrum in accordance with experiment, the theoretical calculations for 
the electron spectrum can be compared with new observations on the intensity of the soft 
component for energies higher than 5 X 10’ ev. The calculations are in satisfactory agreement 


with experimental results found by D. Hall and later by W. Hazen. 


1. INTRODUCTION 


HE problem of the origin of the electrons in 
the lower atmosphere has been considered 
previously by several investigators in the field of 
cosmic rays. The earliest treatments attempted 
' to show that the observed intensity of the cosmic 
rays was a consequence of the interaction of 
high energy primary electrons with the atmos- 
phere, and the total intensities observed were 
secondary products of these high energy pri- 
maries. After the discovery of the mesotron and 
the hypothesis of its disintegration, Euler and 
Heisenberg! assumed that the soft component in 
the lower atmosphere was only partly the result 
of incident primary electrons, and that, in addi- 
tion, there was an appreciable contribution by 
the electrons generated by the disintegration of 
the mesotrons. 

Schein, Jesse, and Wollan,? on the basis of 
experiments carried out in the stratosphere, sug- 
gested that all the soft component might be 
attributed to secondary rays arising from the 
mesotrons—the majority being produced by de- 
cay processes. Later, Rossi and Greisen* calcu- 
lated the intensity of the soft component based 
on this hypothesis and found that the theoretical 
number of electrons arising from the mesotron 
component was in close agreement with observa- 

1H. Euler and W. Heisenberg, Ergeb. d. exact. Natur- 
wiss. 17, 1 (1938). 

2M. Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev. 


59, 615 (1941). 
? B. Rossi and K. Greisen, Phys. Rev. 61, 204 (1942). 
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tions at sea level; but at the altitude of Mt. 
Evans, 14,618 feet, there was an appreciable 
discrepancy between theory and experiment. 
From their results, they concluded that such a 
simple mechanism was unable to account for the 
total soft component in the lower atmosphere. 
However, their method of calculation employed 
integral methods by which only the total number 
of electrons of all energies resulting from the 
mesotron intensity could be found; and the im- 
portant problem of the energy distribution of 
the electrons could not be solved. The approach 
of Rossi and Greisen is also complicated by the 
fact that in calculating the total number of 
electrons, one has to include all the electron 
energies down to very low values where the 
shower theory is less easily treated and becomes 
more complex than at higher energies.‘ Since it 
is considerably more difficult to measure the 
number of electrons with very low energies in 
the soft component, an exact comparison between 
theory and experiment is difficult to achieve. 
Considering the importance of the conclusions 
of Rossi and Greisen, it seemed necessary and 
desirable to try a more nearly accurate and 
direct method for calculating the soft component 
—particularly for the altitude of Mt. Evans 
where the deviations in their calculations from 
the experimental values were about the largest. 
D. Hall has measured the approximate energy 
spectrum of both mesotrons and electrons at this 


* B. Rossi and S. J. Klapman, Phys. Rev. 61, 414 (1942). 
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location, and a more nearly complete theoretical 
solution of this problem can be compared with 
his results. This is possible if the energy spectrum 
of the mesotrons is known immediately above the 
point of observation. A power law spectrum can 
be assumed in accordance with the data on the 
zenith angle distribution and measurements of 
the penetration of mesotrons through lead thick- 
nesses up to 91 cm,} and the calculations which 
follow show that the exponent of the power law 
spectrum can be varied between rather wide 
limits without appreciable change in the spec- 
trum of higher energy electrons. 

The mesotron spectrum was cut off at the 
value of the energy which is necessary to pene- 
trate 18 cm of lead (3.8 10* ev) since the total 
number of mesotrons with energies greater than 
this limit has been measured in balloon experi- 
ments and is known throughout the atmosphere.’ 
Also, below this energy, the power law rapidly 
fails to hold. The mesotrons with lower energies, 
in order to contribute to the soft component, 
must disintegrate within a few radiation units of 
the point of observation; and the electrons 
emitted in these processes will have an energy of 
1.9 108 ev on the average, and will not multiply 
appreciably. Consequently, of these mesotrons, 
those which decay close to the point of observa- 
tion contribute one electron each; and the total 
contribution of the low energy part of the 
mesotron spectrum can be calculated from the 
theory of the disintegration process. 

The theory presented here is based. on the 
cross sections for pair production and brehms- 
strahlung calculated on the assumption of com- 
plete screening, and includes also the ionization 
loss of the shower electrons. This should, then, 
give the correct spectrum of electrons arising 
from the disintegration of the mesotrons for 
energies above 5 X 10’ ev. However, a comparison 
with the experimental values must be considered 
to be only preliminary since the data available 
so far have limited accuracy. The mean life of 
the mesotron used in the calculations has been 
assumed to be the same throughout the entire 


5 It can also be pointed out that the theoretical calcula- 
tions of the mesotron spectrum, based on assumptions 
similar to those used by Euler and Heisenberg, lead to 
this type of spectrum as a close approximation. 


atmosphere, and the value used is 2.15 X 10~* sec. 
as given by Nereson and Rossi.*® 

With suitable changes in the mesotron spec- 
trum, the calculations could be extended without 


modification to give the spectrum of electrons at 


othef altitudes. 
2. CALCULATIONS OF THE ELECTRON SPECTRUM 


In calculating the soft component that arises 
from the mesotron component it is necessary to 
include at least three processes: (a) the decay of 
a mesotron giving rise to a high energy electron, 
(b) the cascade multiplication of electrons, and 
(c) electrons arising from close mesotron col- 
lisions. In addition to these processes, there are 
others (Compton effect and the angular scatter- 
ing of the electrons in the shower) which may 
be considered as corrections although their inclu- 
sion in any other manner makes the mathe- 
matical analysis much more complicated. It can 
be shown, however, that for the lower atmos- 
phere, these corrections are negligible for energies 
above about 5 X10’ ev. 


(a) Decay 


It is assumed that the elementary process of 
the decay of a mesotron consists of spontaneous 
disintegration into an electron and an unknown 
neutral particle (neutrino). It is further assumed 
that the rest mass of the neutral particle is 
negligibly small so that the electron will receive 
on the average one-half of the total mesotron 
energy. According to this hypothesis, in a réfer- 
ence frame in which the mesotron is at rest, the 
electron will be emitted with equal probability 
in any direction and with an energy equal to 
half the rest energy of the mesotron. If the 
mesotron is moving relative to a reference system 
fixed with the earth, it is necessary to use a 
Lorentz transformation to calculate the decay 
in this system. Under these conditions, the energy 
of the decay electron with respect to the earth’s 
system will depend on its angle of emission in 
the system in which the mesotron is at rest, and 
hence on the corresponding angle of emission in 
the earth’s system. If 6 represents the angle in 
the earth’s system between the trajectory of the 
decay electron of total energy E and the trajec- 


* N. Nereson and B. Rossi, Phys. Rev. 64, 199 (1943). 
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tory of the initiating mesotron of total energy E’ 
and rest energy yu, then using the equations of the 
Lorentz transformation one obtains the following 
result: 
2EE’— 
2E(E"?—y*)! 


In Fig. 1, the electron energy E is plotted for 
different values of @ in a polar diagram. The 
polar distance is proportional to the logarithm 
of the total energy E of the decay electron in 
the earth’s system. The various curves corre- 
spond to different mesotron energies. 

Since 108 ev is the critical energy of electrons 
in a shower in air, an electron with this energy 
cannot multiply appreciably in the atmosphere 
and will be brought to rest in a few hundred 
meters of normal air. As a result, these lower 
energy electrons will not contribute appreciably 
to the total intensity of the soft component. 
From the curves in Fig. 1, it is clear that decay 
electrons having energies greater than 3 x 10° ev 
always lie in a small solid angle (less than 11°) 
around the trajectory of the mesotron, and hence 
the approximation introduced by considering the 
decay process as one dimensional is satisfactory. 

The probability that a mesotron will decay in 
a system in which it is at rest is equal to dr/ro, 
where + is the time in that system and 7» the 
mean life of the mesotron. Transformed to the 
earth’s system, this probability beconfes 


cos @ . (1) 


P\(E’)dy=——dy, 


where P;(E’)dy is the probability a mesotron 
with energy E’ will decay in the length dy of its 
trajectory. If the depth from the “‘top” of the 
atmosphere is measured in radiation units ¢ in a 
vertical direction, it follows that 


La p(t) 


thdt=— (2) 


Toe 
where p(t)=a/t represents the barometric for- 
mula expressed in radiation units, and a has the 
value 7.0X10° cm to 7.5X10° cm. This implies 
that the barometric formula in the atmosphere 
is approximated by an isothermal relation. 
Expression (2) gives the number of decay 
electrons generated by a mesotron of energy EF’ 
in a thickness dt at the depth ¢, without regard 
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to its energy. Among these electrons there is a 
probability P2(E’, E)dE that a decay electron 
will have the energy E to E+dE, and on the 
basis of the theory of decay assumed above, this 
probability is given by’ 


PE’, E\dE= (3) 


The probability, then, that a mesotron will 
produce a decay electron of energy E in the 
thickness dt in air is 


P(E’, t)dEdt=P,(E’, t)P2(E’, E)dEdt. (4) 


Fic. 1. Energy of a decay electron as a function of the 
angle of emission. The mesotron trajectory -is directed 
along the line @=0. The mesotron energy E£’ is indicated on 
each curve. For a particular curve, the energy E of the 
decay electron in the earth's system depends on the angle 
@ between the trajectories of the mesotron and electron. 
The length of the radius vector for a given @ is constructed 
proportional to the logarithm of E/10° of the decay electron. 


(b) Cascade Multiplication of Decay Electrons 


The cross sections for brehmsstrahlung and 
pair production used in the process of cascade 
multiplication are taken from the calculations of 
Bethe and Heitler.* These cross sections depend 
in general on the energy of the initiating particle 
(electron or photon). If ¢(£, W)dW is the proba- 
bility that an electron with energy £ will produce 
a photon with energy between W and W+dW, 
and ¥(W, E)dE the corresponding probability for 
pair production, the following expressions apply 
for the cross sections in the limit for high energies: 


brehmsstrahlung, 
Jdv/v; (5) 

pair production, 
¥(u)du =[u?+(1—u)?+u(1—u) (2/3 — 2b) ; (6) 


7 See, for example, reference 3. 
H.A. Betheand W. Heitler, Proc. Roy. Soc. 46, 83 (1934). 
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where E is the electron energy, W the photon 
energy, and v= W/E, u=E/W, with 6=0.0135.° 
These expressions are valid for complete screen- 
ing, and since they do not depend explicitly on 
the energy of the initiating particle, they are 
more tractable mathematically. 

The ionization loss of the shower electrons will 
be included in the calculations. According to the 
Bloch formula, the ionization loss of a particle 
increases very slowly for energies higher than the 
rest energy, and consequently it can be con- 
sidered as a constant a, which has the value of 
10* ev per radiation unit in air. 

At low energies, the Compton effect becomes a 
relatively important factor, but nevertheless this 
process can be neglected in these calculations 
because the energy range considered includes 
only energies higher than 5 X10’ ev. 


(c) Collision Electrons 


Electrons arising from close collision with 
mesotrons will also give a contribution to the 
observed number of electrons and, therefore, 
should not be neglected although the majority 
of the electrons resulting from this source is in 
the lower energy range (below 10*® ev). The 
number of collision or ‘‘knock-on’’ electrons 
in equilibrium with mesotrons has been calcu- 
lated with good approximation by Tamm and 
Belenky.’® Since the average energy of the 
mesotrons at Mt. Evans is much lower than at sea 
level, their calculations show that the contribu- 
tion to the soft component from this process is 
negligible for the electron energiesconsidered here. 

The following calculations were carried out for 
electron energies greater than 5X10’ ev for 
altitudes of 14,618 feet corresponding to a depth 
of 14.3 radiation units below the top of the 
atmosphere. 

Diffusion Equations 

In order to derive the fundamental equations 
for the problem described above, let x(E, t)dE 
represent the number of electrons whose energy 
lies between E and E+dE at an atmospheric 


depth ¢ below a certain depth ft where, sup- 
posedly, — most of the mesotrons are formed. 


ages Roe Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
1T. Tamm and S. Belenky, J. Phys. Acad. Sci. U.S.S.R. 
1, 177 (1939). 
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Similarly y(W, t)dW and f(E’, t)dE’ will be used 
for the corresponding spectra of photons and 
mesotrons, respectively. The principle of con- 
tinuity applied to these differential spectra for 
the present problem leads to the following 
diffusion equation," 


Xo(v)dv+2 o(- 


dr(E, t) 
+a +F(E, t), (7) 
dy(W,t) dv 
=f ), (8) 
where 


F(E, t)= f “P(E', (9) 


1 
D= (10) 
0 
with ¢ measured vertically, and all energies being 
total energies. 


Equations (7) and (8) are the same as the 
diffusion equations of shower theory, except for 
the additional term F(E, #) in (7) which depends 
on the spectrum of the mesotrons. The equations 
can be extended to include the Compton effect ; 
however, its inclusion into the diffusion equations 
complicates considerably the process of solving 
the equations. 

The usual method of solving the diffusion 
equations of the shower theory shows that (7) 
and (8) have unique solutions since the solutions 
not including the term F(£,#) vanish if there 
are no incident particles present at zero thick- 
ness. Also, in case the rest energy of the mesotron 
can be neglected in comparison with the energy 
E’, one can prove that half the total energy 
given out by the mesotrons which disintegrate 
spontaneously appears in shower particles (elec- 
trons and photons).” 

1 The derivation is analogous to that found in H. j: 
Bhabha and W. Heitler, Proc. Roy. Soc. 159, 432 (1937 
J. T. Carlson and J. R. O heimer, Phys. Rev. 51, 220 

H. Phys. 53, 960 (1938); R. Serber, 

hys. Rev. 54, 317 (1938) and summarized in reference 9 

This is only true under the condition that we neglect 


the term representing the ionization loss of the electrons in 
the shower equations. 
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Solutions 


The shower equations (7) and (8) can be 
solved by using Mellin transforms. The homo- 
geneous equations [omitting the term F(E, t)] 
have been studied extensively in the shower 
theory, and the results of these theoretical calcu- 
lations can be used in part in this problem. 

Equations (7) and (8) transformed by the 
Mellin integrals are 


0M,(s, t) 
ot 


=—A (s)M,(s, t) 


+B(s)M,(s, )+Mr(s, t), (11) 
0M,/(s, t) 
———-= C(s) M,(s, t) -DM,(s, t), 


. where the ionization loss a has been assumed 
zero and M,, M,, and My are the Mellin trans- 


(E, ¢t’) 1 +D Jpi(s, E/a) 


[A2(s) +D ]po(s, E/a) 
C(s)g2(s, E/a) 


a 
(s) 
x 0+ f Mp(s, ‘dt — 


M,({s, 0) 


where the functions \;(s) and \2(s) are the roots 
of the equation 


(13) 


the functions p,(s, E/a) and g,(s, E/a) are given 
by 


1 I'(—r)T(s+r+1) 
| pis, f T(s-+1) 
XK r)x~dr, 
| (j=1, 2) (14) 
x r)x*dr ; 


C(s) 
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forms of t), y(W,#), and F(E, respec- 
tively, taken with respect to the energies E 
and W. At the thickness t=0, the transforms are 
denoted by M,(s, 0) and M,(s, 0), and represent 
the boundary conditions. If F(E,t) vanishes, 
(7), (8), and (11) represent the equations of 
showers produced by an incident spectrum of 
electrons or photons. However, if F(E, t) does 
not vanish, the solutions of (11) are similarly 
obtainable by direct mathematical means. It is 
then possible to apply the inverse transformation 
to these solutions which gives the shower func- 
tions (E,¢#) and y(W,#) when the ionization 
loss is neglected. Having found these preliminary 
solutions, it is possible to construct the solutions 
for Eqs. (7) and (8) which include ionization loss, 

With the same independent variable E for the 
electron and photon energies, the complete solu- 
tions, in matrix notation, are™ 


B(s)pi(s, E/a) 
—[A2(s) +D ]gi(s, E/a) 


B(s)p2(s, E/a) ) 

—[Ai(s) +D ]g2(s, E/a) 

| J Mr(s, | 

M,(s, 0) 


and K ,(s, r) is determined by the equation 

A;(s) +D 

XKj(s,r)=rK(s,r—1) (15) 
with the condition 

Ks, 0)=1. 

The paths of integration are lines parallel to the 
imaginary axis, the path of integration c lying to 
the right of the origin and the path c’ to the left, 
both lying in a region where the respective 


integrals converge. 
The formal solution (12) becomes the solution 


+A (s+r) 


43 The notation here is the same as that usually used in 
shower 


theory, see reference 9. 
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" ysually found for showers originating from single 
particles when My =O. If the functions p; and g;, 
which are introduced because of the ionization 
loss of the electrons, are placed equal to unity, 
the solutions go over into those in which the 
jonization loss has been neglected. 

The boundary conditions for the problem con- 
sidered here are 


M,(s, 0)=M,(s, 0) =0, (16) 


which mean that the contribution to the soft 
component from primary electrons and photons 
incident at the top of the atmosphere is negligible 
at a depth of 14.3 radiation units or more. By 
introducing a known spectrum of mesotrons into 
. (12), the line integrals can be evaluated approxi- 
mately by the saddle point method. In all the 
following calculations, only the function (EZ, t) 
which refers to the electron spectrum will be 
considered. 

For the differential spectrum of mesotrons, a 
power law of the following form is assumed, 


f(E’, \dE’=N(t)E’-"dE’, E'>E"” (17) 
with 2.5<n<3.0 and E”=3.8X10* ev.“ The 
function N(t) can be calculated from the meas- 
ured spectrum. Such a law applies closely for 
depths corresponding to Mt. Evans and would 
apply as well or better for smaller depths. For 
greater depths, this simple law will no longer be 
true because a maximum in the mesotron spec- 
trum appears at energies higher than 3.8 X 108 ev, 
this being especially pronounced close to sea 
level. 

As pointed out above, the contribution of 
electrons arising from mesotrons with energies 
lower than 3.8X10* ev to the soft component 
will be small since they will not have sufficient 
energy to multiply appreciably in air, thus 
allowing the spectrum to be broken at this value. 

Introducing the boundary conditions and meso- 
tron spectrum mentioned above into Eqs. (12), 
one obtains for the differential spectrum of the 
electrons the following expression, 


“Such a spectrum is in good agreement with the experi- 
mental data on the zenith angle distribution of mesotrons 
and also with balloon and aeroplane experiments on the 
mesotron intensity. 


N(t’)dt’ 
©) 


rec EE" 2xi Jy 
Au(s) +D 
f (s)+ 
e —Aa(s)] 
X pi(s, ds, 


where the path of integration is a line parallel 
to the imaginary axis. The intercept of this line 
with the axis of reals is such that 0<R(s) <n. 
A similar expression can be obtained for the 
differential photon spectrum. 

Since the mathematical procedure depends on 
the relative values of E and E”, two possible 
cases must be considered. 


(18) 


(a). E<E” 


In this case, the integrand becomes positively 
infinite for s=0 and for s=n, and the method 
of the saddle point can be used. The saddle point 
for the integrand in (18) is determined by the 
equation 


+log (E”/E) 
2s—(n—1) 
(s+1)(m—s) 
which gives the value of s. Using this method in 
evaluating the inner integral in (18), one finds 


finally that the differential electron spectrum is 
given by the expression, 


r(E, t)= 


1 f Nit’) 


2m toc (t' +t) 
Au(s) +D 
(s+1)(m—s)[DAa(s) —A2(s) J 


2st 1)(n—s)+[2s—(n— 
(s+1)*(n—s)?. 
Xpils, E/a)(E" 


(20) 


(0). E>E” 


For these values of E, the method of residues 
can be used since there is a single pole in the 
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20 a0 60 7 
PRESSURE (cm Hg) 


Fic. 2. The mesotron intensity as a function of atmos- 
a pressure, determined by Schein, Jesse, and Wollan. 

dinates represent the number of mesotrons per minute 
received in a vertical direction and able to penetrate 18 cm 
of lead. Abscissae represent atmospheric pressure. 


positive half-plane at the point s=n, and the 
integrand vanishes for large real values of s. 
Consequently, 


pa 
th) =— 
roc 


Ai(n)+D 
(m+1)[Ai(m) —d2(n) 


tN 
x f dt’. (21) 
0 +to 


The expressions (20) and (21) for the electron 
spectrum are considered as the final solutions of 
the problem treated here. Both integrals can be 
evaluated numerically and the corresponding 
spectrum of the electrons found. The actual 
comparison of the magnitudes of the two ex- 
pressions using E”’ = 3.8 X 108 ev shows that (21) 
contributes only a very small part to the integral 
spectrum for values of EZ less than 10° ev. 


Comparison with Experiments 


In Fig. 2 the intensity of the mesotrons is 
shown for an absorber thickness of 18 cm of lead 
as a function of atmospheric pressure.? The 
minimum energy of the mesotrons recorded in 
these experiments was 3.8X10*® ev. With these 
data, the integration of expression (17) will 
give the coefficient of the energy spectrum M(t). 

To calculate the electron spectrum for the 
altitude of Mt. Evans, the numerical quantities 


E/a) 
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TABLE I. Data used in the computations. 


to=2.15X10-* sec. 
u=10% ev 

a=7.2X10° cm 

t= 14.3 radiation units 
n=2.5, 2.6, and 3.0 


E” =3.8X 108 ev 
Vertical incidence 


Mean life of the mesotron 

Mass of the mesotron 

Barometric constant 

Depth 

Exponent (mesotron spectrum) 

Cut-off energy (mesotron 
spectrum) 


listed in Table I were used in the solutions (20) 
and (21). For vertical rays, the differential 
spectrum of electrons as calculated from (20) is 
shown in Fig. 3 for two different values of the 
exponent m in the mesotron spectrum, and for 
three different mesotron cut-off energies E”’. In 
the two lower curves for the energies E”’ = 10° ey 
and E’’=2X 10° ev the value of m was taken as 
2.6, and in the two upper curves, » was taken as 
2.5 and 3.0 with a cut-off energy of E’”=3.8 
X 108 ev. 

The upper curves show clearly that the spec- 
trum of the electrons arising from the decay of 
mesotrons does not depend very strongly on the 
value of the exponent m in the mesotron spec- 
trum. Comparing the areas under the curves, 
one can see that the majority of the electrons in 
the lower atmosphere arise from mesotrons with 
energies lower than 10° ev. The electron spectrum 
rises sharply for lower energies, i.e., most of the 
decay electrons have energies below 10% ev, and 
hence one cannot expect any large number of 


5 
07 


2n 
ENERGY (ev) 

Fic. 3. The differential electron spectrum for an atmos- 
pheric depth of 14 radiation units. Ordinates are x(E, #) 
and the abscissae are the electronic energy E. The curves 
show the changes in character resulting from variation in 
the exponent (m) of the power law for the mesotron spec- 
trum and the cut-off energy (E£”’). 
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ENERGY SPECTRUM OF ELECTRONS 


high energy electrons to be present in the lower 
atmosphere. 

Since these calculations give only the number 
of electrons arising from the disintegration of 
mesotrons of energy greater than 3.8108 ev, it 
isnecessary to include further additional processes 
which contribute to the soft component in the 
lower atmosphere. 

(a). First one must consider the electrons origi- 
nating from the decay of mesotrons with energies 
below 3.8 X 108 ev. To calculate their number, the 
mesotron spectrum for energies below this limit 
must be known. Unpublished measurements of 
Hall give the approximate shape of the mesotron 
spectrum at Mt. Evans for energies below 10° 
ev.» A similar shape of the spectrum is indicated 
in the aeroplane experiments of Schein, Wollan, 
and Grotzinger'® which were carried out at an 
altitude of about 7000 m. Using these data, the 
contribution of decay electrons arising from 
mesotrons of lower energies was calculated by 
means of Eq. (4) taking ionization loss into 
consideration. 

(b). A further contribution to the soft com- 
ponent arises from collision electrons in equi- 
librium with mesotrons. The spectrum of these 
electrons was calculated by Tamm and Belenky.® 
Their data could be used for estimating the 
maximum contribution of knock-on electrons to 
the calculated number of electrons arising from 
the decay. 

The integral spectrum of electrons resulting 
from all the different processes which were con- 
sidered here are represented in Fig. 4. Curve A 
represents the integral spectrum calculated from 
Eqs. (21) and (22) for »=2.5 whereas curve B 
represents the same spectrum for n= 3.0. How- 
ever, the influence of the knock-on electrons in 
the integral spectrum of electrons down to an 
energy of about 5.010’ is negligible as the 
average energy of the mesotrons at Mt. Evans 
is considerably lower than that at sea level and 
practically all knock-on electrons have energies 
below this value. Since experimental data are 
obtainable for energies only down to 7X10’ ev, 
the knock-on electrons can be neglected. This is 


* The author is greatly indebted to Dr. Hall for com- 
munication of his data. 

1%©M. Schein, E. O. Wollan, and G. Grotzinger, Phys. 
Rev. 58, 1026 (1940). 
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FRACTION OF MESOTRON COMPONENT 


3 
ENERGY (x 10” ev) 


Fic. 4. The integral electron s 


trum for an atmospheric 
a of 14 radiation units and vertical incidence. Points 
in 


icated by O represent the results of Hall on Mt. Evans, 
and § represent normalized results obtained by Hazen at 
an altitude of 10,000 ft. 


in very good agreement with recent experiments 
on knock-on electrons by Hazen on Mt. Evans."’ 

If all the electrons in the lower atmosphere 
arise from the mesotron component, then the 
curves of Fig. 4 should represent the measured 
integral spectrum of electrons on Mt. Evans for 
vertical incidence. D. Hall carried out careful 
investigations of the penetration and shower 
production of electrons in different thicknesses 
of lead on Mt. Evans, and his data can be used 
for comparison with the present calculations. He 
determined the number of electrons which pene- 
trated a given thickness of lead whereby the 
electron was identified by its property of pro- 
ducing a cascade shower in about 3 radiation 
units of lead. To obtain the integral spectrum 
from such experiments, it is necessary to calcu- 
late the minimum energy an electron must have 
to penetrate a given thickness of lead. This 
calculation was carried out for three different 
thicknesses (1.5 cm, 3.5 cm, 5.5 cm of lead) 
used in Hall’s experiments.’* By using this 


procedure, Hall’s data for electrons are repre- 


sented by the circles in Fig. 4. All of the ordinates 


17 W. Hazen, Phys. Rev. 64, 7 (1943). ‘ 

18 The minimum energy of an electron which is necessary 
so that at least one particle will emerge below the lead 
thickness was determined by shower theory including 
ionization loss. The values found were: 


Lead thickness Min. energy 
4X 107 ev 


2x 10° 
5.5X 108. 
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56 HENRY E. 
are fractions of the total intensity of the meso- 
trons of all energies at Mt. Evans, which has 
been adjusted to unity. If we use this normaliza- 
tion, the measurements of Hall give the value of 
0.66 for the intensity of the mesotrons passing 
through 18 cm of lead or having an energy 
greater than 3.8X 108 ev. 

As shown in Fig. 4, the theoretical results are 
in approximate agreement with the experimental 
data for the higher electron energies. However, 
below 10° ev the calculations become more com- 
plicated and their reliability is also in question. 
Similarly the experimental evidence for electrons 
with low energies is uncertain because electron 
identification becomes more difficult because of 
the fact that their ability to produce cascade 
showers is greatly reduced. 

Recently Hazen’® has made cloud-chamber 
measurements of the soft component of cosmic 
rays at an altitude of 10,000 ft. By the method 
used, the number of electrons could be estimated 
by counting the showers produced in the eight 
lead plates placed in the chamber. Photons could 
be identified by showers but with the absence of 
an initiating ionizing particle, and mesotrons 
were counted by means of their characteristic 
sharp tracks and the absence of accompanying 
showers. This way, Hazen was able to determine 
the number of electrons with energies greater 
than 108 ev, and from a study of the subsequent 
showers, the energy of each individual electron 
could be determined. 

The results of the experiment showed that the 
total soft component for energies greater than 
2108 ev was seven percent of the penetrating 


19 W. Hazen, Phys. Rev. 65, 67 (1944). 
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radiation. Of this soft component, 2.5/3.5 was 
electronic in nature. Consequently, the total 
number of electrons greater than 2X10® ev is 
about five percent of the total penetrating 
radiation. This value is in good agreement with 
the results of Hall’s experiments and the value 
given by the calculations. Hazen was able to 
determine the energy spectrum of these electrons 
also, and a further experimental check of the 
theory is therefore possible although the com- 
parison can only be considered as approximate. 
The results are shown in Fig. 4 together with 
those found by Hall. 

There are two things which must be con- 
sidered in Hazen’s results. First, the solid angle 
open to electrons in the measurements of the 
soft component was presumably large and since 
zenith angle variations in intensity have not 
been calculated in the present paper, corrections 
for this effect have not been made. Second, the 
difference in altitude between the point where 
the calculations apply and the one used by 
Hazen, which amounts to two radiation units in 
atmospheric depth, is going to modify the meso- 
tron spectrum slightly. In first approximation, 
however, the same exponent in the power law 
of the mesotron spectrum will apply to both 
locations. This means that the electrons gener- 
ated by the mesotrons would have nearly the 
same energy distribution in both cases, their 
total numbers being in direct proportion to the 
total numbers of the mesotrons which produce 
them. 

The writer wishes to express his gratitude and 
indebtedness to Dr. M. Schein for careful col- 
laboration and guidance during the course of this 
investigation. 
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It has been found that when an iron tube is magnetized along its length in a uniform field 
the intensities inside as well as outside the tube are reduced by some amounts which depend 
upon the intensity of magnetization. This is also true for the field outside an iron rod. From 
the distribution of the intensities of magnetization along the tube, which are experimentally 
determined, the amounts of reduction of the magnetic fields may be calculated by applying 
Wiirschmidt’s theory of the demagnetizing factor. The theoretical values are in good agreement 


with those obtained through experiments. 


1. INTRODUCTION 


ALL! found that when a long steel tube is 

placed in a uniform magnetic field the 

intensity of the field H inside the tube is reduced. 
His results are as follows: 

“(1) The reduction of the value of H inside 
the tube becomes greater as the thickness of the 
wall increases. 

(2) The reduction of H within the tube is a 
maximum when the tube is magnetized to a 
condition which about corresponds to the knee 
of the magnetization curve. 

“(3) When the tube is magnetized to a condi- 
tion of saturation, the value of H within the 
tube is practically unaffected by the material of 
the tube.’ 

From these facts Wall concluded that ‘‘one 
aspect of the classical theory of magnetization 
is incorrect,’’® and that ‘‘the magnetization of a 
material is initiated at the surface and proceeds 
inwards by a kind of ‘chain effect,’ so that the 
flux density near the surface is a maximum and 
becomes progressively less as the distance from 
the surface increases.””? 

It was thought that the diminution of H inside 
a tube might be accounted for by the demag- 
netizing effect of the tube itself. Since, unless the 
tube is infinitely long, the applied magnetizing 
field would be affected by the free poles induced 
on its end surfaces, these free poles would have 
different effects at different points in space. 
Consequently, the magnetizing field, and hence 
the intensity of magnetization of the tube, 


1 Wall, Nature 142, 1039 (1938). 
? Wall, Nature 143, 331 (1939). 


would be non-uniform along its length. In 
general, the intensity is high and uniform in the 
central portion and becomes gradually lower and 
less uniform towards the two ends. The range of 
uniform region depends upon the intensity of 
magnetization, being of a minimum at about the 
knee of the magnetization curve and of a maxi- 
mum at saturation. 

The non-uniformity of magnetization appears 
to indicate that the free poles reside not only on 
the end surfaces, but are also distributed along 
the whole length of the tube. All these free 
poles give rise to their own field. The resultant 
field in the neighborhood of the tube is, in 
general, the vectorial sum of three parts, namely, 
the applied field, the field due to the free poles 
induced on the two end surfaces, and that due to 
the free poles distributed along the whole length 
of the tube. The last two are usually collectively 
called the demagnetizing field on account of their 
direction being opposite to that of the applied 
one. For a long tube the effect of the second part 
of the demagnetizing field on the central portion 
of the tube is larger than that of the first. 
Apparently the demagnetizing effect due to the 
distributed free poles other than those on the 
two end surfaces was not taken into consideration 
in Wall’s experiment, although he mentioned 
that the demagnetizing effect was negligibly 
small. Thus it is quite reasonable to assume that 
the reduction of the applied field is simply 
another expression of the demagnetizing effect 
of the whole tube. If this is true, the applied 
field outside the tube should also be reduced. 
It will be seen later that the experimental results 
confirm this prediction. 
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2. THEORY 


According to Wiirschmidt’s* theory, the de- 
magnetizing field of a tube of finite dimension 
can be calculated in the following way. The 
distribution function of the intensity of mag- 
netization of the tube magnetized along its axis 
in a uniform field is given by the relation, 


T=Io{1—C2(x/a)?— C4(x/a)* 


(1) 


where I is the intensity of magnetization at a 
position whose distance from the middle of the 
tube is x; Jo, the intensity of magnetization at 
the middle; a, the half-length of the tube, and 
C2, Ca, Co, etc., the constants to be determined 
from experiments. The free poles induced in the 
section dx at position x from the middle of the 
tube is 


= o( 2x + 


a? 


a‘ 
in which xb? is the cross-sectional area of the 
tube. The free poles induced on the two end 
surfaces are -). As- 
suming the free poles to be concentrated along 
the axis of the tube, we find that the demag- 
netizing field in the central portion, near the 
inner (or outer) surface and in the direction of 


the axis, is 
Cxx*dx 
2Cyx‘dx 


3Cex*dx 
+ 
a?(x?-+ R?)! 


+ f onal 
a*(x?+R?)! 


@ 


where R is the inner (or outer) radius of the tube. 


3J. Wiirschmidt, Theorie des Entmagnetisierungsfaktors 
der Scherun *e von Magnetisierungskurven (Vieweg, 
1925). 


. LIN 


Consider the general term of Eq. (2): 
Du f 
(x?+R?)! 
2n—1 


x 


= (2n— 1) (Dont 2-2), 


(n=1, 2, 3, etc.) 


(x?+ R?)! 


Collecting all the terms for different m’s in suc- 
cession, introducing the limits of integration, and 
neglecting the terms containing R, which is small 
compared with a, we have 


x*dx 
(x?-+ R?)! 


x‘dx 


6x*dx 1 


From Eqs. (2) and (3)—(5), we get 


2xb? 2 83 
Ha= 
a? B 2 


1 
Sot ove 


Hence the demagnetizing field Hz is a function 
of the cross-sectional area, the length of the 
tube Jo, and the constants C2, C4, C6, etc. 


3. EXPERIMENTS 


(a) Determination of the Magnetic Field 
Inside the Tube 


A uniform field was produced by a solenoid of 
120 cm in length and 5.3 cm in inner diameter. 
The iron tube used was an ordinary water pipe 
reduced to a uniform thickness and having a 
length of 60.8 cm, outside diameter of 2.48 cm, 
and inside diameter of 2.29 cm. It was placed co- 
axially with, and in the middle of, the solenoid. 
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IRON TUBE AND ITS MAGNETIZING FIELD 


Fic. 1. The magnetic field H 
inside the tube as a function of 
the exciting current J. 


H in oersteds 


———-—without tube 
with tube 


A search coil was wound on a fiber frame large 
enough to be readily slipped into the middle of 
the tube. The field was measured by the ballistic 
method as usual. 

Before taking readings the tube was demag- 


netized. The exciting current ranging from 0 to 
3.6 amperes was applied in 33 steps. For each 
step the reversal throw of the galvanometer 
deflection was observed, and the mean of the 
two deflections was used in calculation. 


sor 


Fic. 2. The magnetic field H 
outside the tube as a function of 
the exciting current J. 


H in oersteds 


LS 2.0 
I in amperes 


The iron tube was then removed, but the 
search coil kept in the original position. These 
observations were repeated. The magnetic fields 
thus determined are plotted in Fig. 1 against the 
exciting current. 

It is seen that the two curves are similar to 
those obtained by Wall. The only difference is 
that the two curves do not run as close together 
at high exciting currents as do Wall’s. This is 
probably due to the smaller dimensional ratio 


—— — — without tube 
with tube 
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TABLE Ia. Magnetic field outside the rod. : 
— of the tube investigated in the 


(b) Determination of the Magnetic Field 
Outside the Tube 


In this experiment the same solenoid and the 
same iron tube as in (a) were used. The search 
coil is, however, quite different from the previous 
one. It consists of a differential coil with two 
separate windings of equal number of turns 
wound one upon another in opposite sense. Any 
change of magnetic flux inside the inner layer of 
the winding will induce no e.m.f. in the two 
windings as a whole, but any change of flux in 
between the inner and outer layers of the coil 
TABLE Ib. Magnetic field outside the rod. Search coil at Will induce an e.m.f. By inserting a long solenoid 

25 cm from one end of the rod. of a small cross section in the differential coil 

n" and by adjusting the number of turns of the 

Magnetic field in oersteds 4 coil windings, compensation can be rendered 
(without Be complete. 

The field in the middle and on the outside of 
the tube was measured with this differential 
search coil. Measurements were repeated with 
the coil kept in the same position while the tube 
was removed. The results are shown in Fig. 2, 
Curves in Figs. 1 and 2 are remarkably similar. 


(c) Determination of the Magnetic Field 
Outside an Iron Rod 


Since the magnetic field outside the tube is 
reduced, it is to be naturally expected that the 
field outside an iron rod should also be reduced. 
TABLE Ic. M tic field outside the rod. Search coil at A long Armco iron rod, 80 cm in length and 0.6 

13 cm from one end of the rod. cm in diameter, was tested in the same way as 
in experiment (b). Three experiments were per- 
formed. In the first the differential coil was 
placed in the middle of the rod, in the second 
and the third, at 25 cm and 15 cm from one end 
of the rod, respectively. The intensities of the 
magnetic fields thus determined and the amounts 
of percent reduction calculated from them are 
given in Tables Ia-c and the curves plotted 
against exciting current and against exciting 
field are shown in Figs. 3 and 4. 


(d) Determination of the Intensities of 
Magnetization Along the Tube 


For this experiment the same solenoid and iron 
tube as used in experiments (a) and (b) were 
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TABLE II. Distribution of intensity of 


etization along the tube. x = distance from the middle of the tube. a = half- 


length of the tube. J)>=intensity of magnetization at the middle of the tube. J = intensity of magnetization at distance x 
from the middle of the tube. H = magnetizing field in oersteds. 


3 
a 


s/o H=095 1.81 


x/a H=S1.57 64.35 


110.98 119.51 136.62 


o ooo 
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1.00 
1.00 
1.00 
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again used. A search coil of small length was 
prepared, which can slide freely on the tube. 
After demagnetizing the tube, the exciting field 
was increased from 0 to 138 oersteds in sixteen 
steps. For each step the field was kept constant 
and the intensities of magnetization at different 
positions, at 1 or 2 cm apart, were measured. 
Sixteen sets of data are given in Table II where 
the values I/Io are expressed in terms of galva- 
nometer deflections when the search coil was 
placed in the middle and at any other position 
x/a. Eight out of sixteen sets of data are plotted 
in Fig. 5 showing the distribution of intensity of 
magnetization along the tube. 


4. CALCULATION OF THE DEMAGNETIZING FIELD 


The demagnetizing field is calculated from 
Eq. (6) in which the constants c’s are determined 
by Eq. (1). As the terms involving the power of 
x/a greater than 8 are negligibly small, it will 
be sufficiently accurate to use four terms in 
Eq. (1) for calculating the value of J/Jo. For 
the sake of simplicity, replacing any three of 
the values of x/a in Eq. (1) by A, B, and D and 
the corresponding values of 1—I/Io by E, F, 
and G, respectively, we obtain the following 
three simultaneous linear equations: 


Ce, 
BCs, 

Solving for C’s in terms of A, B, D, E, F, and G 


and 


and substituting the actual values of x/a and 
the corresponding values of I/Iy given in Table 
II, we obtain the numerical values of C’s, which 
are given in Table III. The value of Io in Eq. (6) 
is computed on the assumption that the effective 
field within the material is the same as that in 
the free space inside the tube. 

The demagnetizing fields both inside and out- 
side the tube and the corresponding percent re- 


TABLE III. Demagnetizing field inside and outside the 
tube. H;=exciting field inside the tube in oersteds. Ho=ex- 
citing field outside the tube in oersteds. C2, Ci, Cs=con- 
stants. Ha = demagnetizing field inside the tube in oersteds. 
Hao=demagnetizing field outside the tube in oersteds. 
Ha/H;X100= percent reduction of the exciting field in- 
side the tube. Hao/HoX100=percent reduction of the 
exciting field outside the tube. 


(Hai/ Hi), (Heo/ Ho) 
Ha X100 x<100 


duction for different exciting fields are calculated 
and given in Table III. For comparison the per- 
cent reduction of the exciting field obtained 
through experiments and that by calculation are 
plotted in Figs. 6 and 7. It is evident that in 


I/Ie . 
0 1.00 1.00 1.00 1.00 
le 0.07 0.07 1.00 1.00 1.00 1.00 
014 0.14 1.00 1,00 _ 1.00 1.00 
‘h 0.22 0.22 1.00 1 1.00 
0.29 0.29 0.99 1 1.00 
1s 0:36 0.36 0.98 0 1.00 
0.43 0.43 0.97 0 1.00 
Oo 051 0.51 0.96 0 0.99 
058 0.58 0.93 0 0.98 
1S 0.65 0.65 0.89 0 0.96 
0.72 0.72 0.83 0 0.93 
y 076 0.76 0.79 0 0.91 
0:79 0.79 0.74 0.87 
of 0.83 0.52 0.58 0.83 0.69 0.82 0.84 
0.87 046 052 0.87 0.62 4 0.76 0.78 
co) 0.90 0.39 0.44 0.90 0.53 0.66 0.68 
0.92 0.92 0.50 0.63 
in 0.94 031 035 0.94 0.43 0.54 0.54 0.56 
il 095 032 0.29 0.26 0.24 0.24 0.29 0.33 0.95 0.40 O51 0.52 
d 
il 
d 
of 
al 
h 
e 
| Hi Ho Cs Ce Aa 
7 0.95 0.95 058 -0.50 0.81 0.10 0.10 1040 10.17 
181 181 082 -0.79 0.89 0.33 0.32 18.12 17.65 
353 3.52 108 1.35 1.13 1.10 32.09 31.19 
6.11 6.10 109 1.23 259 251 42.43 41.15 
8.70 8.70 1.09 1.37 1.32 3.98 3.85 45.67 44.27 
17.37 17.21 0.88 0.97 6.51 6.18 3746 35.91 
25.92 25.70 0.69 -O.52 0.76 7.03 6.84 27.14 26.64 
38.79 38.55 0.37 0.09 0.38 6.02 5.91 15.52 15.32 
$1.62 51.36 0.23 0.09 0.50 5.20 5.12 10.08 9.98 
64.50 64.20 0.08 0.27 O41 4.11 408 637 6.36 
73.03 72.76 0.06 0.16 0.54 4.10 3.87 S61 5.32 
81.78 81.44 —0.004 0.24 0.46 3.34 334 408 4.10 
93.51 9453 0.004 0.10 0.57 3.25 3.26 348 3.44 
110.24 111.71 0.020 —0.08 0.69 3.20 3.19 2.90 2.86 
118.69 120.32 0.001 0.06 0.66 3.28 3.28 282 2.72 
135.49 137.75 0,004 0.68 2.96 2.96 2.19 2.45 


H.in oersteds 


lin amperes 


Fic. 3. The magnetic field H outside the rod as a function of the excitin 
current J. a, the search coil is in the middle of the rod; 6 and c, the searc 
coils are at 25 cm and 15 cm from one end of the rod, respectively. The 
upper J axis is for curve a, the middle for curve b, and the lower for curve c. 

he dotted line represents the field without the rod; the solid line represents 
the field with the rod. 


each case the theoretical and experimental curves hardly be explained by Wall’s theory of mag- 
nearly coincide, except in the central part where _ netization. It is merely a necessary consequence 


the former is a little higher than the latter. of the theory of demagnetization; for the free 
poles should produce their demagnetizing effect 
+ See in the whole space. The striking similarity of the 


The reduction of the magnetic field outside two pairs of curves, shown in Figs. 1 and 2, 
the tube or rod as shown in Fig. 2 or Fig. 3 can strongly support this argument. Furthermore, it 


Per cent reduction of H 


0 20 40 60 80 100 10 
Hin oersteds 


Fic. 4. The nt reduction of the exciting field H outside the rod as a func- 
tion of the exciting field H. A, the search coil is in the middle of the rod; B and C, 
the search coils are at 25 cm and 15 cm from one end of the rod, respectively. 
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IRON TUBE AND ITS MAGNETIZING FIELD 


Exciting Field 
Oereveds 
0.95 


4 


Fic. 5. The distribution of intensity of magnetization along the tube. 
I is the intensity of magnetization at x cm from the middle of the tube; 
Ie Se ay of magnetization in the middle of the tube; a, the half-length 

the tube. 


is of interest to compare the behavior of the about 20 cm long, is entirely negligible; (2) the 
curves shown in Fig. 3. It will be noticed that demagnetizing effect due to the free poles dis- 
those curves in solid lines coincide in the case of _ tributed along the rod is also negligibly small in 
a and 6 with those in broken lines at the ends of _ the central region at high exciting currents. We 
the curves. This means: (1) the demagnetizing are thus enabled to understand that the magnetic 
effect of the end surfaces on the central region field within the tube is practically unaffected by 


Per cent reduction of H 


0 30 60 90 120 
H in oersteds 


Fic. 6. The percent reduction of the exciting field 7 inside the tube as a 
function of the exciting field H. The broken line gives the experimental 
curve and the solid line gives the theoretical curve. 
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Fic. 7. The percent reduction of the exciting field H outside the tube 


as a function of the exciting field H. The broken line gives the experimental 


curve and the solid line gives the theoretical curve. 


the material of the tube when magnetized to 
saturation as claimed by Wall. If the exciting 
current does not exceed 2.2 amperes, the demag- 
netizing effect on the central region produced by 
the free poles cannot be overlooked. The effect 
reaches a maximum when the exciting current is 
within the range of 0.2 to 0.5 ampere. On Wall’s 
theory it seems difficult to explain why the per- 
cent reduction of the magnetic field represented 
by curve C should be higher than that by 
curve B, and that by curve B should be higher 
than that by curve A (Fig. 4), but these can be 
easily understood from the demagnetizing effect 
due to the free poles distributed along the whole 
length of the tube, which falls to a minimum in 
the middle of the rod and increases towards its 
two ends. 

An inspection of the general features of the 
family of curves shown in Fig. 5 will render the 
above statement clear. Two points are worthy 
of note: (1) The slope of each curve increases, in 
general, with the value of x/a, indicating that 


' the density of free poles and hence the demag- 


netizing fields produced by them also increases 
from the middle towards the end of the tube. 
Although the free poles may be few in number 
in the central portion of the tube, they may have 


notable demagnetizing effect in their immediate 
neighborhood. Unless the tube is sufficiently 
long, this effect cannot be neglected. (2) When 
the exciting field increases from zero onward, 
the slope of the curves for the same value x/a, 


but not too close to x/a=1, also increases and © 


reaches a maximum value, corresponding to the 
exciting field of 8.7 oersteds; any further increase 


of the exciting field makes the slope of the curves — 


gradually decrease. There is reason to infer that 
if the tube is long enough and the exciting field 
high enough, the slope of the greater part of the 
curve will decrease to zero as its limit, i.e., the 
magnetization of the greater part of the tube is 
uniform. Since the slope of the curve at a given 
point denotes the density of the free poles in 
the corresponding position of the tube, and the 
density of the free poles is proportional to the 
demagnetizing force, we may say that the de- 
magnetizing force of the free poles increases from 
zero to a maximum value corresponding to the 
exciting field of 8.7 oersteds and then decreases 
when the exciting field increases from zero up to 
a considerably high value. The above statement 
qualitatively explains, on the theory of demag- 
netization, the nature of the curves shown in 


Figs. 1-3. 
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6. CONCLUSION 


From the above discussion and from the close 
agreement between the calculated percent re- 
duction of the magnetic field and the experi- 
mental values shown in Figs. 6 and 7, it may be 
concluded that the reduction of the magnetic 
field inside and outside the tube is due to the 
demagnetizing effect of the free poles distributed 
along the whole length of the tube, or in other 


words, the reduction of the magnetic field in the 
neighborhood of the tube can be accounted for 
by the theory of demagnetization without intro- 
ducing any other hypothesis. 


ACKNOWLEDGMENT 


The author wishes to thank Dr. J. W. Shih 
for many helpful discussions in the course of this 
investigation. 


PHYSICAL REVIEW 


VOLUME 66, NUMBERS 3 AND 4 


AUGUST 1 AND 15, 1944 


On the Existence of Atomic Nitrogen in the Upper Atmosphere of the Earth 


Ta-You Wu 
Department of Physics, National University of Peking, and Institute of Astronomy, 
Academia Sinica, Kunming, China 
(Received June 10, 1944) 


Evidence for the existence of atomic nitrogen in the upper atmosphere is presented. It is 
shown that the absence of certain lines of NI and the absence of continuous absorption which 
leads to the dissociation of Nz cannot be considered conclusive evidence against the existence 
of atomic nitrogen. Relative densities of the E and F layers are discussed, and it is shown that 
observations of high density in the F layer may be compatible with the existence of atomic 
nitrogen. Possible transitions which produce the night sky spectrum are discussed. 


ARTLY on the basis of the observed spectra 
of the night sky and of the aurora and partly 
on the basis of the absence in the laboratory of 
continuous absorption by the nitrogen molecule 
leading to dissociation N2+hy-N-+N, the view 
is held by many investigators that in the upper 
atmosphere oxygen exists in the atomic and 
nitrogen in the molecular state. On this view 
many attempts have been made by various in- 
vestigators to deal with such problems as the 
existence of various ionized layers, the density 
at different heights in the atmosphere, the excita- 
tion mechanism of the radiations of the aurora 
and the night sky, etc. The purpose of the 
present note is to present some arguments for the 
view that there may exist atomic nitrogen in 
the upper atmosphere, and to indicate the 
changes that this picture entails in our present 
view of the above-mentioned problems. 


One of the reasons for supposing nitrogen to 
exist in the molecular state in the atmosphere is 


that while the forbidden lines 1D—1'S, *P—'D of 
oxygen OI are the most prominent radiations in 
night sky and the aurora, the forbidden lines 
*S—?P,4*S—?D of NI are absent in these spectra. 
Kaplan has pointed out the possible presence of 
the ‘S—?P line in the spectrum of the night sky." 
While this identification, if confirmed by more 
observational data, would establish the existence 
of atomic nitrogen, it must be emphasized that 
the absence of this and the ‘S—?D line constitutes 
no proof against the presence of atomic nitrogen. 
The transition probability of ‘S—*P is of the 
order of 1 per second and the lifetime of the *P 
state is thus comparatively short. Any N(?P) 
atoms produced by dissociation by solar radia- 
tion (see Section II below) will rapidly settle 
down into the states *?D and *S on the emission 
of the lines *7D—*P and *S—?P. The lifetime of 
the *D state is of the order 10° seconds so that 
the failure to observe the line ‘S—*D may only 
mean an insufficiently large number of N(?D) 
atoms. Unless N(?P) are excited by proper and 


1 J. Kaplan, Phys. Rev. 56, 858 (1939) ; 57,'249 (1940). 
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frequent processes in the night sky and in the 
aurora, there may be considerable N(‘S) and 
N(D) atoms without revealing themselves in 


the spectrum. 
II 


Perhaps the main reason for believing in the 
absence of atomic nitrogen in the upper atmos- 
phere is the absence of continuous absorption by 
N2 down to 600A in the ultraviolet that leads 
to the dissociation of the molecule. Now, accord- 
ing to the configuration theory, of the many 
theoretically possible molecular states resulting 
from two N atoms, the following can combine by 
dipole transitions with the normal state X ',+ 
of N2: two states from N(?7D)+N(?D); one 
and three 'II, states from N(?D)+N(?P); 
and two '2,*+ and three 'II, states from N(?P) 
+N(P). Some, if not all, of these are repulsive 
states so that continuous absorption by Nz in 
the normal state X 'Z,+ will lead to dissociation. 


Thus 

N2(X My), (1) 
N2(X (?P)+N(2D) (42+; Mu), (2) 
N2(X 12,+) +hyN (2D) +N (2D) (I). (3) 


That such continuous absorptions have not been 
observed in the laboratory may be due to the 
circumstance that the potential energy curves of 
these repulsive states rise very steeply with the 
decrease in internuclear distance so that, by the 
Franck-Condon principle, absorption will not be 
' appreciable until the radiation frequency is very 
much greater than vo given by hv»p= D+E, where 
D is the dissociation energy of the normal state 
X '%,+ and is 7.35 volts, and E is the sum of the 
excitation energies of the two resulting N atoms 
from the ‘4S state. Thus appreciable absorptions 
by the processes (1), (2), (3) may lie in the region 
beyond 500A not yet explored in the laboratory, 
and in the region \>600A the absorption may be 
so weak that they are masked by the very strong 
absorptions of the Hopfield and the Worley- 
Jenkins bands that lead to the ionization of N2. 

Despite the situation that dissociative absorp- 
tion will take place only for radiations of very 
short wave-length, an appreciable number of N2 
may still be dissociated since there is now strong 
evidence that there is a great abundance of high 
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frequency radiations emitted by the sun.? They 
are probably due to the radiative recombinations 
of electrons and ionized atoms carried through 
the photosphere by convective transfer, and 
their energy density may be very much greater 
than that of a blackbody at ~6000°K, possibly 
by some high powers of 10. As shown later, a cal- 
culation of the degree of dissociation of Ne is yet 
not possible on account of the lack of pertinent 
informations. The point, however, is that the 
processes (1), (2), (3) are theoretically possible 
provided there are radiations of sufficiently high 
frequency, for the existence of which there seems 
now to be strong astrophysical evidence. 


Consider now briefly some of the consequences 
of the existence of atomic nitrogen in the upper 
atmosphere. On account of the smaller mass of 
the N atom, the dissociation of Nz will lead toa 
greater density, at a given height, than that 
obtained on the assumption of a Ne and O 
mixture. Thus estimates based on the assumption 
that Nz and O are in diffusive equilibrium above 
~150 km show that the air density in the F 
layer (~250 km) is about 1/5000 of that in the 
E layer.’ However, from the observed collision 
frequencies by radio wave experiments, it is 
known that the air density in the F layer is 
about 1/30 of that in the E layer. This ab- 
normally high density in the F layer may be 
partly due to the dissociation of Nz into N atoms. 

The dissociation of N: will also give rise to an 
ionized layer due to the ionization of the N 
atoms. When the variation of the degree of 
dissociation with height is known, it will be 
possible to calculate the height of this ionized 
layer. This may contribute to a better under- 
standing of the nature of the various ionized 
layers now known. Unfortunately a calculation 
of the variation of the degree of dissociation with 
height is not yet possible for reasons given later. 


IV 


Another observed fact connected with the 
problem of the ionized layers may perhaps be 


um M.N.R. A.S. 101, 177 (1941); 102, 
3S. K. Mitra and H. Rakshit, Ind. J. Phys. 12, 47.(1938); 
S. K. Mitra and A. K. Banerjee, ibid. 13, 117 (1939). 
‘Bates, Buckingham, Massey, and Unwin, Proc. Roy. 
Soc. A170, 322 (1938). 
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understood on the introduction of atomic nitro- 
gen in the upper atmosphere. The effective re- 
combination coefficients of electrons are observed 
to be ~~1.210-* and 0.8 cm/sec. in the 
E and F layers, respectively. In the prevalent 
theory,® the electrons are removed by direct 
radiative recombinations with positive ions 


nt++e—n+hyv, recombination coeff. =a,, 


and through attachment by a neutral particle 
forming a negative ion, and subsequent non- 


- radiative recombinations between the positive 


and the negative ions, yielding two excited atoms 


n+e—n-+hv, attachment coeff.=a, 
n~+nt—n'+n", transfer coeff. = a;. 


On denoting the concentrations by simply e, n-, 
nt, and letting 
h=n~/e, 


then, since e+n-=n*, one has 
(1+A) (de/dt) =dn*/dt, 
dn+/dt= 
de/dt= — (ae+)a,)e’. 


Theoretical calculations of Bates et al. show that 
for oxygen, cm/sec. and a,~10-" cm?/ 
sec.> Thus to obtain an effective recombination 
coefficient in agreement with the observed values 
above, the ratio \ must be of the order 100- 
1000 in the E layer and 10 in the F layer. The 
high concentration of negative ions in the E 
layer is supported by the observed variations of 
terrestrial magnetism and is understandable if 
one assumes that there oxygen atoms are abun- 
dant and, having an electron affinity of 2.2 
volts, readily form O- ions. The equilibrium 
ratio \ of negative ion — electron is determined by 
the rates of attachment above and of detachment 


O-+X—-OX-+e, detachment coeff.=d, 
where X is a neutral particle, namely, 
O- a0 


and 


The summation is taken over all the probable 
detachment processes, for example, X =O, Oz, No, 


‘Bates, M. N. R. A. S. 100, 25 (1939), has calculated 
the theoretical absorption coefficient of N atom. 


etc. Now on the prevalent view of a Ne, O 
atmosphere, the concentration of O would in- 
crease relative to that of N2 in passing from the 
E to the F layer so that a decrease of \ in 
passing from the E to the F layer is difficult to 
understand. On the other hand, on assuming Nz 
to be dissociated, the ratio O/N of the concentra- 
tions of oxygen and nitrogen atoms will decrease 
with height. Since the N atoms are probably 
just as active in the detachment process but form 
no negative ions on account of their very small 
electron affinity, the value A will accordingly 
decrease in passing from the E to the F layer. 

Bates et al. have pointed out the possibility of 
thermal detachment in the F layer as a cause of 
the comparatively low concentration of negative 
ions. As the electron affinity of O atom is 2.2 
volts, the temperature must be quite high for 
thermal detachment to be important. As the Nz 
are mostly dissociated by radiations of frequency 
very much greater than » given by hyx=D+E 
(see Section II above), the resulting atoms will 
possess high kinetic energies. This would result 
in a high temperature in the upper atmosphere, 
on which there seems to be general agreement 
among different authors but which has been 
ascribed to the energy of electrons in the ioniza- 
tion of the atoms rather than of the atoms in 
the dissociation of the molecules. 


Vv 


A difficulty of the prevalent view that nitrogen 
exists as N» and is hence confined to the lower 
strata (below ~250 km) of the atmosphere is 
that, according to it, the observed nitrogen bands 
from high aurora become difficult to understand. 
There does not seem to be a plausible theory 
that can account for the excitation of the various 
radiations of nitrogen, especially for the sunlit 
aurora that have been observed at heights as 
great as ~800 km and of which the negative 
bands of N.* are particularly enhanced. The 
theory of Maris and Hulburt attempts to ascribe 
the aurora to particles sent up with great veloci- 
ties to very great heights on collisions of the 
second kind with excited atoms or molecules in 
the uppermost region of the atmosphere, ionized 
by solar radiations, and spiraling down along 
the earth’s magnetic lines of force toward the 
polar regions. It has been shown, however, 


recently by the writer that the mechanism sug- 
gested by these authors is entirely inadequate for 
explaining the aurora.* On our present view here, 
the extension of atomic nitrogen to great heights 
may at least provide for the presence of nitrogen 
molecules in the upper atmosphere by such 
recombination processes as 


N+N—-N2+hy 
(Vegard-Kaplan and other bands of N:), 


(first negative bands of N-*). 


VI 


Finally, on the assumption that there exist in 
the upper atmosphere atomic oxygen and nitro- 
gen produced by the processes (1), (2), (3) and 


O2(X (4) 


the writer has recently proposed a theory of the 
excitation mechanism of the night sky spectrum.’ 
The key processes are the radiative recombina- 
tions between the atoms, forming metastable 
molecules of and O:: 


N(*S) +N ?2,*; *2u+)+hy, (5) 
*2u*)+hy, (6) 
*Zut)+hy. (7) 


The strong Vegard-Kaplan bands are excited by 
(5) and (6). Collisions among the various meta- 
stable atoms and molecules provide for the 
excitation of the various radiations. For example, 
the green and red lines of OI are excited by the 
processes 


20,(C *2)0.(X *2)-++0(¢P) 
+0O(!S)+0.12 volt, (8) 
O(?P) +N (2D)—O(!D) +N (4S) +0.47 volt. (9) 


A high metastability of O2(C *2) will enable its 
concentration to build up and the process (8) 
will explain the small rise in intensity of the 
green line '!D—'S in the first part of the night. 
On the other hand, with an initially large number 
of O(@P) and N(?D) atoms produced by (1), (2), 
(3), (4), and their gradual depletion by recom- 


6 T. Y. Wu, Proc. Ind. Acad. Sci. A18, 345 (1943). 
7T. Y. Wu, Proc. Ind. Acad. Sci. A18, 40 (1943). 
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bination and other processes, the process (9) 
may explain the decrease of intensity of the red 
line *P—1D immediately after sunset. Without 
assuming the presence of N atoms, the excitation 
of the Vegard-Kaplan bands and the other 
bands of Nz seem very difficult to explain.’ 


Vil 


From the foregoing discussions, it seems very 
desirable to settle the question whether N, is 
dissociated in the upper atmosphere by calcu- 


lating theoretically the degree of dissociation as — 


a function of the pressure in the atmosphere. 
This would necessitate a knowledge of the 
radiation energy density p, of the proper wave- 
length, the absorption coefficient 7, as a function 
of frequency, the temperature 7 of the atmos- 
phere where dissociation is being calculated, and 
the effect of ‘‘cycles’” of such processes for the 
formation of the normal N2(X 'Z,+) and the 
removal of N(4S) as 


*2,+) 
N2(A 12,*) 
+Vegard-Kaplan bands. 


On ignoring the effect of such processes and con- 
sidering the process (3) alone, one readily obtains 
the following dissociation formula 


 2ghs (MAT? 
e-DIkT (| — 
[N2] 8x Sah 


f / f Tye y2dy, 
vo ¥0 


where [N ], [N2], and the g’s are the concentra- 
tions and the statistical weights of the atoms and 
the molecule; w is the vibrational frequency, 1 
the nuclear separation, 2M the mass of the 
molecule N2(X '!2); W is the dilution factor and 
is ~4X10-*; and p is the partial pressure of N; 
at the height where the degree of dissociation is 
being calculated. As [N_?/[N2] is exceedingly 
sensitive to the values of p,, JT, and 7,, and as 
these quantities are not known with any cer- 
tainty, any estimate of the degree of dissociation 
with the use of questionable values for these 
quantities will be premature. 
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The quantum theory of radiation describes the effect of retardation in the interaction of 
two electrons by the mutual emission and absorption of light quanta. A study has been made 
of a process involving the interaction of two electrons and any number of quanta. It was noted 
that the number of terms involving virtual quanta needed to describe the retardation in the 

' usual first-order approximation of perturbation theory depends on the number of real quanta 
involved in the process under consideration. It was shown that, nevertheless, the resulting 
expression obtained by summing over all the intermediate states involving virtual quanta is of 
the same form for any process. Mgller’s formulas for the retarded interaction of two electrons 
with no light quanta and with one light quantum emitted are special cases of our equation. 


I. INTRODUCTION 


T is well known that the quantum mechanical 

Hamiltonian of a system of charged particles 

obeying Dirac’s equation and the radiation field 
can be written! 


H=H,+H,+H,.+H', (1a) 


where H, represents the Hamiltonian of the 
particles, H, that of the light quanta, H, the 
electrostatic energy of the particles, and H’ the 
interaction of the particles and the field. We have 


H,= (1/2) > ik (1b) 
H’ = — ex(ax-A(k)), 


where ¢; is the charge and a, the Dirac matrix 
for the kth particle, A(k) the vector potential 
evaluated at the position of e, and ry the dis- 
tance between the ith and th particle. 
Equation (1a) contains an explicit expression 
for the instantaneous Coulomb interaction only, 
as given by (ib). The effect of retardation is 
included in H’ and is described by a mutual 
emission and absorption of light quanta. In the 
simplest possible case of retarded interaction of 
two particles in the absence of radiation, it 
consists of the emission of just one quantum by 
one particle and reabsorption by the other and 


and 


* A thesis submitted in ial fulfillment of the require- 
ments for the degree Doctor of Philosophy in the 
Faculty of Pure Science, Columbia University. 

t Publication assisted by the Ernest Kempton Adams 
Fund for Physical Sesmat of Columbia University. 

‘For the general theory and the notation used see 
W. Heitler, The Quantum Theory of Radiation (Oxford 
University Press, 1936), Chapter III. 
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vice versa.? The purpose of this paper is to show 
that while in the more general case of the inter- 
action of any number of light quanta with two 
particles the retardation corresponds to emissions 
and absorptions described by a large number of 
terms involving virtual quanta, we can sum all 
these terms and obtain an expression for the 
retarded interaction which is of the same form 
for any process. 


Il. THE MATRIX ELEMENT 


As it is not possible with the present methods to 
solve the Schrédinger equation Hy = —th(dy/dt) 
with the Hamiltonian (1a) directly, a perturba- 
tion theory is required. The method often used 
treats the terms (1b) and (1c) for the interaction 
as small perturbations causing transitions from a 
state Ey with a certain energy and momentum 
of the particles and a certain number of. quanta 
to another state Ey. The transition probability 
is zero unless the total energy is conserved, and 
then it is given by the well-known equation*® 


= (24/h)pg,| (2a) 


where pg, is the density of the final or initial 
states with energy Eo, and 


is the matrix element for the transition. The 


? This case was first calculated with a correspondence 
method by C. Moller, Zeits. f. Physik 70, 686 (1931) and 
Ann. d. Physik 14, 531 (1932). The present method was 
first used for this problem by H. Bethe and E. Fermi, 
Zeits. f. Physik 77, 296 (1932). 

* W. Heitler, reference 1, p. 90. 
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matrix element has a value different from zero 
for the Coulomb interaction of two particles, in 
which case it is equal to 


(3a) 


for two free particles, and for transitions involv- 
ing emission or absorption of just one light 
quantum, where it is given by 


eet»). (3b) 


Here the wu’s are the Dirac amplitudes of the 
electrons and a, is the component of @ in the 
direction of polarization of the light quantum ). 
The momentum is conserved for these transitions. 

If the matrix element for the direct transition 
is zero, we have to use a succession of simple 
transitions leading from the initial to the final 
state through a number of intermediate states, 
for which we still have conservation of mo- 
mentum, but no longer conservation of energy. 
We have to replace (2b) by 


H'= 
(Eo — E1)(Eo—E2) (4) 


X (Eo— Er-2)(Eo—Er-1) 


where the denominator contains the product of 
the differences of the energies of the initial and 
intermediate states, and the numerator the 
products of the matrix elements for the suc- 
cessive simple transitions. These transitions have 
to be chosen in such a way that the total ex- 
pression (4) is of the lowest possible order in the 
electronic charge e for which it does not vanish,‘ 
but otherwise we have to sum over all possible 
successions of intermediate states. We can see 
from (3a) and (3b) that the matrix element for 
interaction with a light quantum is proportional 
to e and the Coulomb interaction is proportional 
to e? and therefore equivalent to two interactions 
with light. 

After these general considerations let us now 
consider a process involving z+Z light quanta 
and two particles. Then our energy balance for 


the initial and final states is 


‘It is a well-known difficulty of the quantum theory of 
radiation that any calculation of a process in a higher 
than the lowest possible order leads to divergencies. 
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where £,°, E,” and E,°, E,¥ are the energies of 
the first and second particle in the initial and 
final state, respectively, and a---z and A---Z 
are the energies of the light quanta. We have 
written the equation in a form corresponding to 
absorptions only, but our calculations include 
emissions as well, as the only difference would be 
a change in sign throughout for the energy and 
momentum of the particular light quantum. 
Furthermore, we can include positrons as well as 
electrons by changing the sign of the energy and 
momentum of the particle. 
The conservation of momentum requires 


Pi’ +p2°+at+b+ 
(6) 


The first non-vanishing matrix element for such 
a process is of the order z+Z-+2 in e. It involves 
z+Z transitions due to the interaction of the 
particles with all the light quanta and one 
transition due to the instantaneous Coulomb 
interaction of the two particles. We obtain a 
matrix element of the same order if we have 
instead of the instantaneous Coulomb interaction 
a retardation described by one transition corre- 
sponding to the emission of a light quantum K’ 
by one particle and a second transition due to 
the absorption of this quantum by the second 
particle. Any succession of transitions involving 
the same light quantum more than once or more 
than one Coulomb interaction would either be 
of a higher order than those just considered or 
would not fulfill the conditions (5) and (6). 

We note that we cannot expect (4) to give 
finite results for any values of the energies and 
momenta of the particles and quanta involved 
subject to the conditions (5) and (6) only. We 
assume for the moment that we have to deal 
with non-vanishing denominators of (4) only and 
shall discuss the limitations of the first-order 
approximation later. 

Let us now consider the possible successions 
of interactions in detail. We can have one 
particle interacting with any m quanta, and the 
other one with the remaining z+Z—m quanta. 
Furthermore, the interactions of the m quanta 
with the first particle can take place in any order 
(and likewise for the z+Z—mn quanta and the 
second particle). In addition the Coulomb inter- 
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action can take place after the interaction with 
any »+M quanta. Then one has to sum over 
both signs of spin and energy of the electron in 
the intermediate states. Finally the exclusion 
principle has to be taken into account by making 
the total matrix element antisymmetric in the 
coordinates of the two particles.® 

We do not attempt to sum over all these 
possibilities since we want only to study the form 
of the retardation for a given Coulomb inter- 
action. From now on we limit our calculations to 
that part of the sum, which has the same nu- 
merator in (4), namely, 


KT (ait anus”) a yu”) 
11 n+2, M+2 


(2 Ap|? for the 
instantaneous Coulomb interaction 


7 
X Max /K'] 


for the retardation, 


with K’= | Ap|, the absolute value of the change 
in momentum of either particle due to the 
Coulomb interaction. 

For the moment we consider the two parts for 
instantaneous and retarded interaction sepa- 
rately. The part containing the instantaneous 
Coulomb interaction corresponds to the following 
process : 


Electron e; emits (or absorbs) 
a; b; ooo N 

Electron e2 emits (or absorbs) 
A:B;---M 


Instantaneous 


€, emits (or absorbs) 
n+1; 

emits (or absorbs) 
M+1;-:- Z. 


We denote the light quanta (and their energies) 
by small letters if they interact with the first 
electron and by capital letters for the second 
electron. 
We have already arbitrarily ordered the quanta 
which interact with particle 1, and separately the 
‘For this last point compare, e.g., N. F. Mott and H. 


S. W. Massey, The Theory of Atomic Collisions (Oxford 
University Press, 1933), Chapter XV. 


Coulomb interaction 


quanta which interact with particle 2. We still 
have to sum over all possible orders of inter- 
action subject to these limitations. For example, 
we can have one process where é; interacts with 
a, b, --+ m and then e: with A, B, --- M but also 
another one where e; interacts with a, then e2 
with A, then again e; with }, c, --- and similar 
possibilities. The only restriction (beyond the 
arbitrary partial ordering of quanta mentioned 
above) is that the order of interaction for each 
particle has to be the same and that the Coulomb 
interaction can take place only after e,; has 
interacted with the light quantum m and é, with 
M, but before their interactions with n+1 and 
M-+1, respectively. 

We note that the part of our sum which con- 
tains the instantaneous Coulomb interaction can 
be factored into a product of the Coulomb term 
and of two series, the one containing all the 
successions of transitions taking place before the 
Coulomb interaction and the other containing 
the successions of transitions taking place there- 
after. 

Now we have to consider the case where the 
Coulomb interaction is replaced by the emission 
and subsequent absorption of a virtual quantum 
K’, which transfers the same amount of mo- 
mentum. It is essential to keep in mind that 
this need not take place in two consecutive steps.® 
The only limitation is that e, can only emit K’ 
after its interaction with n and before its interaction 
with n+1, and that e2 can only absorb K’ after its 
interaction with M and before its interaction with 
M-+1 and vice versa. The extreme case would be 
the succession of interactions where one particle 
has emitted K’ and interacted with all z quanta 
and is therefore in its final state, while the 
second particle is still in its initial state. 

We note that we can no longer factor the sum 
describing the retardation as in the case of 
the instantaneous Coulomb interaction. Further- 


' more, we see that because of the fact that the 


emission and absorption of K’ need not take 


* This point has not been emphasized before. The ex- 
pression given by W. Heitler, reference 1, p. 101, Eq. (31), 


takes account of consecutive ste 
the general matrix element for 
two particles. The same criticism applies to the paper by 
W. Heitler and L. Nordheim, J. de phys. et rad. 5, 449 
(1934), which gives the same expression. The actual calcu- 
lations of this paper neglect the retardation, however, so 
that no error is introduced. 
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place in two consecutive steps, the sum de- 
scribing the retardation will contain more terms 
than the sum containing the instantaneous inter- 
action. The number of terms needed to describe 


the retardation depends on the number of quanta. 


involved in the entire process. This appears 
strange, at least formally, as there seems to be 
no reason why the expression for the interaction 
of one particle with radiation and for the instan- 
taneous interaction of two particles should be of 
the same form independent of the process, but 
not the one for the retarded interaction. But we 
have to keep in mind that we are still describing 
the retarded interaction by virtual quanta, and 
we shall see that the expression obtained after 
eliminating them is again independent of the 
process under consideration as it should be. 

We first simplify the notation for the energy 
differences appearing in the denominator of (4). 
If the first transition were due to the absorption 
of light quantum a by é:, then this particle would 
have the energy E;*=(|pi°+a|*+ 7)! and the 
difference between the energy of the initial and 
the first intermediate state would be 


a+E,°—E;* 
as the second particle is unchanged. We denote 
this difference by (E,*). Similarly we put 
(8a) 
(8b) 
Ex" 


If now we have an intermediate state reached 
by the interaction of e; with m quanta and of é2 
with Z quanta, then the energy difference is 


=(E:")+(E2"). (8c) 
Il. A THEOREM 


with 


Let 
Ao=0, Ai, ‘Ay 


Bo=0, Bi, Ba, 


1 By 
BitAnsy: Angi(BitA 
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be two sets of numbers, subject only to the re- 
striction —B, for any 7, 

Let (A, B) stand for any product of n+m 
factors of the form A;+B;, where we obtain 
each factor from the preceding one by replacing 
either A or B by the following number of its set, 


Theorem: The sum of the reciprocals of all] 
possible (A, B) obtained by starting with either 
Aot+B, or Bot+Ai is 


We shall prove the theorem by induction. It is 
certainly true for »=1 and m=1 since 


1 1 1 
Ai(A B,(Ai+B;) A,B, 


As the sum does not depend on the initial order 
of the numbers in the two sets, it is sufficient to 
show that if we add another number A»: at any 
place to the set A,---A, and construct S,.;,,, 
it will be equal to Sy,m/Ans:. The proof can be 
demonstrated most easily if we insert A,,, 
between Ao and A;. We shall use the symbol 
(A: +B,)—(Aisi to indicate a step which 
changes an A to the value indicated, and like- 
wise for a change of the B. 

We now proceed to construct Sy+1,m. We start 
with the factor (Bo +Ani1). Summing all the 
1/(A, B) which can be obtained by (Bp +Ani1)> 
(Bo+A,) and then proceeding according to our 


rule, we simply obtain where sit 
indicates those terms of Sy, which start with 
(Bo+A;). For 
we obtain S,, m/(Bi+An4:1) and for Ant) 
we obtain BiSa,n/ 
[A n41(BitA n41)_] where indicates the sum 
of those terms of S,, m which start at (B,+Ao) and 
continue with A unchanged up to (Ao+B,) and 


_ then change to (A,+8B,). 


Proceeding in this fashion and adding the 
B; 
terms with the same S,,» as factor, we obtain 


A 


B 
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and the general term is from the method of construction 


B,B,---B, 


Bz 


where S;,z indicates the sum of all the (A, B) 
formed from Ao=0, Ans: and Byo=0, By, Bz---B, 
and for which the theorem holds true by assump- 
tion since x<m and 1<n. 

Adding the right-hand side we obtain 


Al Bi Bo Bz 
Sat1,m= Sant: + Sa. m+ eee 
+Sn,m)/A nti=Sn,m/Anti 
which is what was to be proved. 


IV. SUMMATION OF THE SERIES 


We first sum the series defined in Part II 
containing the instantaneous Coulomb interac- 
tion. We remember that we could factor this 
series and we have already defined our (E,”) and 
(E2") by the relations (8) in such a manner that 
we can immediately apply our theorem to the 
first factor containing the successions of inter- 
actions taking place before the instantaneous 
Coulomb interaction, i.e., up to (E:") and (E™). 

Before we can sum the part containing the 
interactions taking place after the instantaneous 
Coulomb interaction, we have to simplify the 
notation for this part also. But we cannot use 
the same definitions as before since now the 
momentum of the particles has been changed not 
only by the interaction with the light quanta, 
but by their Coulomb interaction as well. The 
change in momentum of the first particle is 
equal and opposite to the change in momentum 
of the second one, and is given from (6) by 


Ap=p."*!—p." =p." —[p:°+a+b+---+z] 
= (9) 
We define now 


(10a) 
Ey 

(E:°)= (10b) 
-Q—---—Z|?+ 


where now (£;?) is the energy difference for an 
intermediate state in which the first particle has 


Bz 
San = 


Bz 
Sa m/A n+ls 


interacted with the first p—1 quanta and the 
second particle has interacted with all Z quanta 
and similarly for (E2°). The energy difference for 
an intermediate state in which e,; has interacted 
with p—1 quanta and e: with Q—1 quanta, is 
given by 
a+b---+(p—1)+A+B::- 

which from (5) is equal to 

Now we can apply our theorem again if we let 
(Ey*), (E:**")-- -(E,"*") correspond to 
and (E27-")- -(E,™ +) to B,---By. 

Hence we obtain finally for the sum of the 
terms containing the instantaneous Coulomb 
interaction, omitting for the moment the factors 
common to the entire series given in (7), 


| Ap|? 


1 
(Ev 


We now proceed to the summation of the part 
of the series containing the retardation. We 
noted before that whereas the terms containing 
the instantaneous Coulomb interaction could be 
separated into the product of a series containing 
terms up to m and M only and another series 
containing the terms from +1 and M+1 up, 
no such separation is possible for the terms con- 
taining the retardation. Our former definitions for 
the (E)’s were consistent as we always had to 
deal with (E)’s from the same part of the series. 
Now, however, we may have an intermediate 
state where the first particle has already emitted 
the (m+1)th light quantum while the other 
particle has not yet absorbed the K’ quantum 
(describing the retardation), i.e., we have terms 
of the form 


(11) 
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and if we define now 
AE=+E," 
(12) 
we can write this term 
(13a) 
For the case where K’ had been emitted by the 


second particle, we obtain a general term of the 
form 


which can be written 
+(E2%)-AE—K’. (13b) 


Let us sum first the part corresponding to the 
emission of K’ by the first particle. Calling 
AE—K’'=X, we have then a series where each 
term is the product of expressions defined as 
before, except that at the step (E,")—>(E,"*") we 
do not take (E£,"*'), but (E£,"+")+X, then 
(E,"**)+X and so on. However, the step 
(E2”)—(E.”*) subtracts the X again. The sum 
can therefore only include those terms where 
(E"*") is reached before (E2”*). 

We split the sum into partial sums according 
to the point at which the K’ is absorbed, i.e., 
into sums with the same r occurring in the transi- 
tion where 
r can have any value between +1 and z+1 
(with (EZ,*+')=0). Then each sum can be factored 
into two sums containing the terms before and 
after this transition, respectively. Then we can 
apply our theorem without difficulty up to 
this transition, with (Z,*), ---(E."), +X, 
(E,"**)+X, ---(Ei")+X corresponding to Aj, 
--+A, and (E24), ---(E2”) corresponding to 
- ++ Bm. We can furthermore apply it to the terms 
after this transition, with (E;*), ---(E.") corre- 
sponding to Ai, ---A, and (E27), ---(E;“+") to 
Bi, -+-Bm. (For the single term obtained for 
p=z+1 we have simply (E27), - -(E2”**).) 


Adding the partial sums again we obtain 
1 
-(E1*) 


1 


(Ey"*") +X) - (Ey) tee 


1 
((Ey"*1)+-X) ((Ey***)+-X) eee ((E1*) 


But the expression in brackets is just equal to 
the sum 5°1/(A, B) formed from the sets 0, X 
and 0, (Z;*)---(E,"*), and therefore the entire 
sum is equal to 

We can sum the second part, corresponding to 
emission of K’ by the second particle, in exactly 
the same way, only having —AE—K’ instead of 
AE—K’. Adding these two sums we obtain 


K’ 


2K’ 
Now we have to sum over the directions of 
polarization of K’ and add the expressions (11) 
and (14). As the independent directions of polari- 


zation are mutually perpendicular and perpendicu- 
lar to the direction of propagation of K’, we have 


= . 


Remembering that K’= Ap, we obtain finally 
for that part of the matrix element (4) which 
corresponds to the interaction of z quanta with 
e, and Z quanta with é2 with fixed order of inter- 
action for each particle and which corresponds to 
Coulomb interaction after m and M as discussed 
before 


(14) 


(—e)*+2+2[ 
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where Ap and AE are given by (9) and (12), 
(Ex*)::*(Ei") and (E24)---(E:”) by (8) and 
+ and - -(E3”) by (10). 

Equation (15) includes the summation over all 
the terms of (4) having the same numerator (7), 
but the summation over all the different possible 
numerators has still to be done for any particular 
process as discussed before. 
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V. DISCUSSION 
Comparison of the expression containing the 
retarded interaction (15) with the expression 
containing only the instantaneous Coulomb in- 
teraction (11) shows that the sole difference is 
that the factor 


|Ap|* 
is replaced by 


independent of the process under consideration. 
This is in spite of the fact that in (4) the number 
of terms needed to describe the retardation did 
depend on the process. Equation (16) is simply 
the usual matrix element (3a) for the instan- 
taneous Coulomb interaction applied to the 
particular transition n—n+1, M->M-+1 of our 
calculations. We see that by eliminating the 
virtual quanta, we have obtained an expression 
(17) analogous to the matrix element (3a) which 
can be used as a matrix element for the retarded 
interaction for any process involving the inter- 
action of radiation and two electrons.’ Using (17) 
the retarded interaction can then be treated as a 
single transition in the succession of intermediate 
states appearing in the total matrix element (4). 
Therefore the above noted difficulty of a de- 
pendence of the form of the retarded interaction 


on the particular process under consideration was__. 


only apparent. 

While we have obtained this correspondence 
of the expressions (16) and (17) for the instan- 
taneous and retarded interactions from the 
summed expressions (11) and (15), we can never- 
theless, if desirable, use our ‘‘matrix element of 
retarded interaction’’ even before summing with- 
out introducing an error. The formalism of 
perturbation theory assigns a meaning to the 
summed expression only, and we are therefore 
justified in using any expansion whatsoever, as 
long as the summed expressions agree. But we 

7It should be remembered that our result has been 
derived for a process involving two particles. Its validity 


for processes involving more than two particles remains to 
be proved. See also footnote 6. 


could not have obtained the correspondence 
without first summing our series. If we fix our 
attention on just one term of the series involving 
the instantaneous interaction, i.e., a fixed order 
of interactions of the light quanta with the first 
and second particle, it is impossible to assign a 
correspondence to a definite part of the series 
involving the retardation. This is due to the fact 
that, as pointed out before, the emission and 
absorption of the virtual quanta need not take 
place in two consecutive steps. 

In the case that no radiation is present (15) 
reduces to (17) alone. Equation (17) can then be 
further simplified. Ap is just equal to pi’ 
and we obtain from Dirac’s equation a(p;" —p,°)u 
= —E,°)u=AE-u and therefore our expres- 
sion reduces to 


| Ap|?—(AE)?* 


which is just the expression obtained by Méller.? 

It can be seen similarly that (15) reduces to 
the expression obtained by Moller for the inter- 
action of two particles and one light quantum,’ 
the only other case for which the matrix element 
has been calculated directly. In this case we have 
still to include a summation of the expressions 
(15) over the four possibilities that the retarded 
interaction can take place before or after the 
interaction with the light quantum and that the 
light quantum can interact with either one of the 


§C. Mgller, Proc. Roy. Soc. 152, 481 (1935). 
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two particles. In addition to these cases (15) can 

be used to calculate the emission of several 
quanta at the collision of two particles. It can 
also be applied to the problem of creation of an 
electron-positron pair due to simple or multiple 
scattering of a y-ray, as in the Dirac theory the 
presence of a positron in the final state can be 
described by the presence of an electron with 
negative energy in the initial state. 

It should be noted that while we started from 
an expression which involved both particles in 
intermediate states simultaneously, each factor 
in the denominator of expression (15) involves 
one particle in an intermediate state at a time 
only. Therefore the method of Casimir® for 
summing over both signs of the energy of the 
particles in the intermediate states can now be 
applied directly. 

Furthermore we can easily see that our ex- 
pression (15) reduces to the known formula for 
the case where one particle can be considered as 
infinitely heavy. Then the energy is conserved 
for the interaction of the electron and the quanta, 
but not the momentum, because of the presence 
of the nucleus. The heavy particle cannot inter- 
act with the light quanta, but only with the 
electron. Therefore it does not pass through 
any intermediate states, and we have to put 
and (u2*au?)=0. With these 
changes, (15) together with the definitions (8a) 
and (10a) correctly describes the interaction of 
one electron with radiation or the creation of a 
pair of particles by the scattering of y-rays in 
the field of a nucleus. 

In the course of our calculations we have 
arbitrarily assumed that none of the energy 
differences in the denominator of (4) vanishes, 
and we have now to study the limits of applica- 
bility of (15). This expression diverges in two 
cases. The first case occurs when any of the 
(E;*)’s or (E2”)’s vanish (found by putting the 
expressions (8a, b) and (10a, b) equal to zero). 
As each of these (Z)’s refers to one particle only, 
the limitation introduced hereby is the same as 


*H. Casimir, Helv. Phys. Acta 6, 287 (1933). 


HAVAS 


for a process involving the interaction of radia. 
tion and a single particle. In the familiar cases 
this excludes processes which involve light 
quanta with the same energy and momentum in 
both the initial and final states and which could 
therefore be described by an approximation of 
lower order in e than the one applied. 

The second case occurs when | Ap|*—(AE)*=0, 
This is the case for vanishing interaction of the 
two particles. Then our process splits up into 
two independent processes each involving a 
single particle and a certain number of quanta. 
Again these could have been calculated in a lower 
order approximation. This case is similar to the 


situation for the calculation of the scattering of 


two charged particles in a non-relativistic first 
order approximation, where we obtain the 
Rutherford formula, which diverges for zero 
deviation, i.e., vanishing interaction of the two 
particles. There, however, the exact solution is 
known, and it again reduces to the Rutherford 
formula for all but vanishing deviation angles, 
for which it also remains finite.’ Unfortunately 
the exact solution of the relativistic two-body 
problem is not known. 

Inspecting formulas (8c), (10c), and (13a, b), 
we might expect still other cases of divergencies 
if these expressions vanish, while these cases 
cannot be seen directly from our final Eq. (15). 
But these expressions can be indefinitely small 
without any indication of a maximum in the 
summed expression, due to cancellations. This 
can be seen directly from the proof of our 
theorem. If a number A; of one set is almost 
equal to the negative of a number B, of the other 
set, we can still apply the theorem. If we 
call B,=—A;i+e, the sum depends upon ¢ as 
1/(—A?2+A.«), with e=0, which is a monotoni- 
cally varying function of «. Therefore (15) does 
not show any resonance-like behavior in the 
neighborhood of the values obtained by putting 
the expressions mentioned above equal to zero. 

The author wishes to thank Dr. Joseph M. 
Keller for his interest and helpful discussions. 


10 W. Gordon, Zeits. f. Physik 48, 180 (1928), in par- 
ticular Eq. (7’), p. 190, and the following discussion. 


PH 


= 

me 

Ell 

res 

fou 

2.5 

sec 

anc 

in é 

anc 

228 

aga 

1.9 

1.9 

of I 

val 

net 

full 

tha 

of 1 

whi 

the 

Fun 

2) 

3, 

‘ 


PHYSICAL REVIEW 


VOLUME 66, NUMBERS 3 AND 4 


Lifetimes of Resonance Lines of Cadmium* 


Haro_p W. WEBB AND HELEN A. MESSENGER** 
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The mean lifetime of the triplet resonance line 5'So—5*P1, 3261A, of cadmium was measured 
by the alternating voltage method. The value found was 2.14 X 10-* second, with an estimated 
precision of three percent. Similar measurements made to determine the mean lifetime of the 
singlet resonance line, 5'S9—5'P:, 2288A, showed that the excitation which resulted in the 
emission of this line was mainly due to cascading from a higher level (or levels) by infra-red 
transitions, the excitation to the 5'P; level by electron impact being relatively much less 
probable than that to the higher level. The lifetime of the infra-red transition was found to be 
9.0 X 10-* second. From the form of the curve showing the dependence of the persistence of the 
radiation emitted by the excitation tube upon the vapor concentration, the mean lifetime of 
the 2288A transition was estimated as approximately 2.1 10~-* second. The lifetime of the 
subordinate series lines 5086-4800-4678A was found to be 10~7 second, but since the measure- 
ments were affected by absorption and re-emission of the radiation, this value can be taken 
only as the maximum possible value of this lifetime. For the diffuse subordinate series lines, 
53P—5D, similar measurements showed that their lifetime was not greater than 10-7 second. 
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- lifetimes of each of the two resonance 
states of cadmium, 5'P; and 5*P;, has been 
measured by several methods. Koenig and 
Ellett,! by measuring directly the decay of the 
resonance radiation in a stream of cadmium vapor, 
found for the lifetime of 3261A (5'S)—5*P), 
2.5X10-* sec. Kuhn? obtained the value 2.3 x 10-° 
sec. by the magneto-rotation method, and Ellett* 
and Soleillet,* from studies of the depolarization 
in a magnetic field, obtained the values 2.3 X 10-° 
and 2X10-* sec., respectively. The lifetime of 
2288A (5'So—5'P:) as measured by Kuhn,5 
again using the magneto-rotation method, was 
1.98X10-* sec. Zemansky*® found the value 
1.99 10-* sec., using the line absorption method 
of Ladenburg and Reiche. Soleillet* computed the 
value of this lifetime as 10-* sec. from the mag- 
netic depolarization data; however, to account 
fully for his results, he resorted to the assumption 
that there existed a second hyperfine energy level 
of the 5'P, state with the lifetime of 10-* sec., 
which is not in accord with the generally accepted 
theory. 


* Publication assisted b Le the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

** Hunter College, New York, New York. 

1H. D. Koenig and A. Ellett, Phys. Rev. 39, 576 (1932). 

*W. Kuhn, Naturwiss. 13, 725 (1925). 

3A. Ellett, +. Rev. 33, 124 (1929). 

+P. Soleillet, omptes rendus oo 212 (1928). 

ow. Kuhn, Naturwiss. 14, 48 ( 1926). 

*M. Zemansky, Zeits. f. Bhs 72, 587 (1931). 


Measurements of the lifetimes of these two 
resonance Hines of cadmium made by use of the 
alternating voltage method, previously employed 
to measure the lifetimes of lines of hydrogen, 
mercury, and potassium,’ are described in this 
paper, together with some results obtained in 
studying the lifetimes of several subordinate 
series lines in the cadmium spectrum for which 
no earlier measurements are on record. 


METHOD 


The method used in determining the mean 
lifetimes of lines has been described in the several 
earlier papers’ and will be but briefly discussed 
here. The radiation studied was produced in an 
excitation tube which contained cadmium vapor 
at low pressure and which was provided with a 
grid for controlling the excitation. This radiation, 
isolated by suitable filters, passed to a photo- 
electric tube so designed that saturation occurred 
at very low values of its control-grid voltage. An 
alternating sinusoidal voltage was applied to the 
control-grids of the respective tubes in phase. 
The bias voltage on the excitation tube grid in- 
sured that excitation occurred only during a 
known fraction of the positive half-cycle, and a 


7F. G. Slack, Phys. Rev. 28, 1 (1926); H. W. a ont 
H. A. Messen; er, Phys. Rev. 33, 319 (1929); 
Randall, Phys. Rev. 35, 1161 (1930) ; P. H. Garrett, Meare 
Rev. 40, 779 (1932); D. Sinclair and H. W. Webb, Phys. 
Rev. 50, 440 (1936). 
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suitable bias voltage on the photoelectric tube 
caused the current to the collecting electrode to 
reverse with the voltage on the control-grid. At 
low frequencies radiation reached the photo- 
electric tube only in the positive half-cycle, and 
hence there was no reverse photoelectric current. 
At higher frequencies, the radiation, because of 
its persistence, reached the tube in both half- 
cycles; the positive current was reduced, and 
reverse current appeared. The net photoelectric 
current to the collecting electrode therefore 
diminished with increasing frequency and ap- 
proached asymptotically a ‘final level’’ when the 
persistence resulted in a uniform emission from 
the excitation tube. 

The first step in the measurement of r, the life- 
time of a line, by this method, was the determi- 
nation for this radiation of the ‘‘static charac- 
teristics’ of the excitation tube (radiation as a 
function of grid voltage) and of the photoelectric 
tube (photoelectric current as a function of grid 
voltage, for constant radiation). The alternating 
voltages were then applied to the control-grids of 
the two tubes in phase and the photoelectric cur- 
rent measured for a suitable range of values of 
the frequency f. The voltages applied to the grids 
were made the same for all frequencies to main- 
tain constant the total radiation excited. The 
results are expressed as the ratio R of the 
photoelectric current at the frequency f to the 
photoelectric current at a frequency so low that 
the persistence of the radiation was practically 
negligible. This will be called the “reference fre- 
quency” hereafter. Data from the curve obtained 
when R was plotted as a function of f-were used to 
compute 7, the mean lifetime of the line. 

In measuring the lifetimes of resonance lines 
the absorption of the radiation presents a serious 
complication. Each photon absorbed produces an 
excited atom which persists in its excited state 
with the same mean lifetime as an atom originally 
excited by electron impact. The observed radia- 
tion consists, then, of photons which come 
directly from atoms excited by impact and 
photons which leave the excitation tube after one 
or more absorption-re-emission processes. This 
increases the persistence in the mass of vapor in 
the excitation tube, and hence the ‘‘apparent”’ 
lifetime T is greater than 7, the difference in- 
creasing with the vapor pressure. To determine r 
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it is therefore necessary to correct for the effect of 
absorption and re-emission, which it is difficult to 
compute precisely, or else to make measurements 
at pressures so low that absorption is negligible. 
This occurs at such low pressures that the radia- 
tion is usually insufficient for precise measure- 
ments. For example, in measuring the lifetime of 
2288A, it was necessary to work at a vapor pres- 
sure below 1.5X10-* mm. At this pressure the 
kinetic mean free path of the electrons in the 
vapor was of the order of 10* cm, so that the 
efficiency of the excitation tube was exceedingly 
low. This difficulty, usually encountered only in 
measurements with resonance lines, proved 
troublesome also in measurements on subordinate 
series lines when intense excitation was used. 


APPARATUS 


Figure 1 shows schematically the essential 
elements of the apparatus used, for the measure- 


cw 
Oscillator 


Fic. 1. Diagram of apparatus. 


ments with 3261A and 2288A. In the main it was 
similar to that used in earlier experiments and 
fully described by Garrett.* The excitation 
tube is shown in cross section. The envelope was 
of fused quartz; S was an indirectly heated, oxide 
coated cathode; G, the control-grid, was of nickel 
mesh and covered only the front half of the 
cathode; W, the wall grid, was of nickel mesh on 
the front, or window, side and of solid sheet in 
the rear. Since it was essential that absorption of 
radiation be reduced to a minimum, the tube was 


* P. H. Garrett, reference 7, p. 781. 
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so designed that by the application of suitable 
voltages to the grids excitation would occur 
mainly between G and W and radiation reach the 
window by traversing only a thin layer of vapor. 
The tubes were sealed off after evacuation in the 
usual manner, with especial care to eliminate 
traces of mercury. 

Adjoining the excitation tube was a reservoir 
V, on the walls of which a thin layer of carefully 
purified cadmium metal had been evaporated 
before the tube was sealed off from the pumps. 
This reservoir was immersed in an oil bath, 
electrically heated, and maintained by a thermo- 
stat within 0.5°C of the temperature corre- 
sponding to the vapor pressure desired in the 
tube. The body of the tube was in an oven con- 
sisting of a copper cylinder with very heavy walls 
heated by gas jets. A fused quartz window oppo- 
site the cathode transmitted the radiation. The 
temperature of the oven was held only slightly 
higher than that of the oil bath to permit the 
concentration of cadmium atoms WN in the tube 
to be taken as that corresponding to the tempera- 
ture of the bath without correction. This concen- 
tration was computed from the vapor pressure 
given by the formula in the International Critical 
Tables for the appropriate range. 

The photoelectric tube used in studying the 
resonance radiations is shown in cross section in 
the figure. The envelope was of quartz. The 
sensitive surface was magnesium, selected for the 
measurements with 3261A and 2288A because its 
long wave-length limit is in the neighborhood of 
4000A and its maximum sensitivity is near 
2500A.* This was deposited on the plate P, a 
nickel sheet 1 cm X4 cm. The control grid H was 
rectangular in cross section and was of nickel 
mesh in the front and of sheet nickel on the sides 
and back with a clearance of 2 mm between it 
and P. The shape and clearance were designed to 
give current-voltage (H — P) characteristics which 
would reverse sharply with the voltage and to cut 
down the time lag of the tube produced by the 
time of passage of the photoelectrons between P 
and H. The auxiliary electrode K was also of 
nickel mesh and served to prevent the return to 
P of electrons which had passed through the grid 


“SAL. H L. Hughes and L. A. DuBridge, Photoelectric Phe- 
nomena (McGraw. Hill Book Company, Inc., New York 
and London, 1932), p. 162. 


H into the region between H and K. The plate P 
was coated with magnesium by the method de- 
scribed by Rentschler.* An auxiliary electrode 
(not shown in the figure) was used, consisting of 
a magnesium wire, 1 mm in diameter, placed 
between H and K, parallel to the axis of the tube. 
After baking the tube and outgassing the elec- 
trodes, the tube was filled with argon at a suitable 
pressure and with the magnesium wire as cathode 
and the electrodes P and H as anode, a coating of 
magnesium was deposited on P by sputtering. 
The argon was then removed by absorption by | 
cocoanut charcoal contained in a side tube 
immersed in liquid air. After the argon had been 
absorbed, this tube was sealed off. 

Provision was made for interposing between 
the excitation and photoelectric tubes various 
filters including a quartz cell containing pure 
cadmium vapor. A water cell, with windows of 
quartz or of glass with high ultraviolet trans- 
mission, was interposed between the tubes when 
it was necessary to protect the photoelectric tube 
from the heat of the oven. 

Photoelectric current was measured by a sensi- 
tive quadrant electrometer. To avoid disturb- 
ances on opening the electrometer key and to 
simplify the high frequency circuit, the lead from 
P was connected to ground through a condenser 
C of approximately double the capacity between 
P and H. The electrometer was connected in 
parallel with this condenser through a resistance 
R of 13X10° ohms, which isolated it from the 
high frequency circuit. 

A vacuum-tube oscillator with loosely coupled 
tuned secondary circuit was the source of the 
alternating voltage. A calibrated peak voltmeter 
was employed to measure and monitor the 
voltage. By the use of switches, not shown in the 
figure, either grid could be connected inde- 
pendently to sources of direct potential or to the 
source of alternating voltage. When alternating 
voltage was used on both grids, their terminals 
were connected as directly as was possible to 
minimize any differences between the magnitudes 
and phases of the voltages on these two electrodes. 

The cadmium absorption cell employed as a 
filter in some of the measurements was especially 
useful in determining the purity of the resonance 


*H. C. Rentschler and D. E. H Frank. Inst. 
223, 135 (1937). 
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radiations. It was a cylindrical cell of fused 
quartz, 2 cm in diameter with plane faces 1 cm 
apart. Solid cadmium in a side tube, the tempera- 
ture of which was controlled by an electric heater, 
furnished cadmium vapor at the desired concen- 
tration. The body of this cell was enclosed in an 
oven to keep the vapor at a temperature slightly 
higher than that of the side tube. Quartz windows 
permitted the transmission of the radiation to 
and from the cell. 


EXPERIMENTAL PROCEDURE AND RESULTS 
3261A 


The measurements of the mean lifetime of the 
resonance triplet resonance line 3261A were 
made with vapor concentrations in the excitation 
tube ranging from 8.7 X10” atoms cm=* (204°C) 
to 3.310" atoms cm~* (280°C). The radiation 
2288A, which was very intense at these concen- 
trations, was filtered out with a Corex D filter 
(Corning Glass Company), which absorbed all 
wave-lengths shorter than 2500A. When testing 
for the presence of wave-lengths other than 
3261A in this filtered radiation, the cadmium 
absorption cell was interposed between the ex- 
citation tube and the Corex D filter. With G—S 
equal to 11.0 volts, and the excitation tube at 
214°C, approximately the temperature at which 
precise measurements were desired, no measurable 
amount of radiation was transmitted by the 
absorption cell when the vapor pressure therein 
was 0.05 mm (300°C), showing that of wave- 
lengths above 2500A the excitation tube was 
emitting substantially only 3261A. The absorp- 
tion coefficient for this line computed from the 
absorption measurements taken at lower pres- 
sures in the absorption cell, was 1.310-“N. 

Measurements to determine the lifetime were 
made with the bias voltage W—S=4.5 volts, and 
bias voltage K—P=4.5 volts. In each determi- 
nation of the lifetime it was first necessary to 
determine the ‘‘static characteristics” of the ex- 
citation and photoelectric tubes for the radiation 
in question at the selected vapor pressure. When 
the intensity of 3261A radiation in the excitation 
tube was studied as a function of the G—S 
voltage, taken with H—P constant at three volts, 
the curve first showed the excitation of 3261A at 
G—S=4.5 volts, with no correction for contact 
potential or drop in the leads. The accepted 


excitation potential for this line is 3.8 volts. At 
higher voltages there were several breaks in the 
curve, which departed sufficiently from a linear 
relation between radiation and voltage to necessi- 
tate a graphical calculation when determining an 
accurate value of the lifetime. The static charac. 
teristics of the photoelectric tube were of ex- 
ceptionally simple form. Saturation was reached 
when H—P was little more than 0.4 volt. The 
reversal point was at H—P=-—0.1 volt, making 
it unnecessary to use a bias voltage on H—P 
during the alternating voltage measurements, 
The reverse current was approximately 30 percent 
of the positive saturation current and varied only 
5 percent as H—P was varied from —2 to —9 
volts. In such measurements special care is neces- 
sary to eliminate the effects of charges which 
accumulate on the walls of the tube and give 
erroneous results. 

When the alternating voltages were applied to 
the control-grids, the peak value used was 6.6 
volts, with a fixed bias of 4.5 volts applied to 
G—S, so that radiation was excited only in the 
positive half-cycle. The crest voltage attained by 
G—S, 11.1 volts, was sufficiently low to insure 
that the radiation excited and transmitted 
through the filter was substantially 3261A and 
that the ionization was negligible. A test made 
with alternating voltage applied to G—S only, 
H—P being held at constant voltage, showed the 
electrometer current to be independent of the 
frequency, indicating that the voltage applied at 
the external terminals of G and S was the same 
as that between the electrodes G and S and that 
the behavior of the excitation tube could be 
predicted from the static characteristics. With 
the alternating voltage applied directly to G-—S 
and H—P, a series of measurements was made to 
determine the R—f curve for each of several 
concentrations of cadmium vapor in the range 
specified above. In this case R was the ratio of 
the photoelectric current at the frequency f to 
that at the reference frequency, 1300 cycles. It 
was found that the final level R,, the value of R 
approached asymptotically by the curve at high 
frequencies, agreed well with that calculated 
from the static characteristics. This is an im- 
portant check, signifying that the behavior of the 
photoelectric tube could be correctly predicted 
from these characteristics. 
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The lowest concentration at which the energy 
was sufficient to make precise measurements 
possible was that corresponding to 218.5°C. At 
this concentration the secondary radiation re- 
sulting from absorption and re-emission which 
reached the photoelectric tube was less than 10 
percent of that resulting from direct excitation. 
At this concentration the R—f curve was taken 
with special care, and a precise value of R for 
f=80,500 cycles, 0.647, was determined by re- 
peated observations. Then from the R—f curves 
taken at oil bath temperatures of 244°, 232°, 215°, 
and 203°C, respectively, the values of the fre- 
quency fo, for which R had this same value, 0.647, 
were also determined. Graphical calculations 
made from the respective static characteristics of 
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Fic. 2. Curve to determine lifetime of 3261A. r=1.826xf 
X10-* sec. for zero concentration of cadmium atoms. 


the excitation and photoelectric tubes show that 
if the radiation had been emitted from the ex- 
citation tube according to a simple exponential 
relation (i.e., all radiation direct and no absorp- 
tion and re-emission), 7, the lifetime of the line, 
would have been given by r= (1.8262f9) sec. In 
Fig. 2 the values of 1.826xfo, determined from the 
five curves, are plotted as a function of the vapor 
concentration. Extrapolation of this curve to 
zero concentration gives the correct value of 1/r 
since the effect of absorption is thus eliminated. 
This extrapolation would ordinarily be difficult to 
carry out with precision because the form of the 
curve for low concentrations is not easily calcu- 
lated from the tube geometry. However, for the 
tube used, the form of the curve is approximately 
known from the very careful measurements of the 
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lifetime of the 2537A line of mercury made by 
Garrett,* who used an excitation tube approach- 
ing closely in size, shape, and positions of elec- 
trodes the tube used in these measurements. The 
curve drawn in Fig. 2 is the curve of Fig. 3b of his 
paper, with the scales so adjusted that the curve 
passes through the precisely determined point for 
218.5°C and best approximates the experimental 
curve for higher values of the concentration. The 
departure of the experimental points from this 
curve is in every case within the experimental 
error of the measurements. The ratio of the scales 
of abscissas is approximately the reciprocal of the 
ratio of the atomic absorption coefficients of the 
two lines. For zero concentration the curve gives 
1/7 =0.467 X 10° sec... The mean life- 
time of the line 3261A is therefore determined as 
2.14X10-* sec., the precision of the measure- 
ments being estimated as better than three 
percent. 

This result is about seven percent lower than 
values given by the magneto-rotation and mag- 
netic depolarization methods.2** The value 
2.5X10-* sec., found in the measurements of 
Koenig and Ellett,! was probably subject to 
larger error than the other methods because of 
the difficulties inherent in the method. 

In the measurements here reported, the excita- 
tion of the 5*P, level by cascading from higher 
levels was considered as a possible source of error. 
However, a series of studies made with varied 
peak voltages indicated this excitation to be 
negligible. 

Some additional measurements made with 
3261A are of interest because of the light they 
throw upon the effect of the absorption-re-emis- 
sion process. These were made with a second 
excitation tube which was practically identical in 
construction with the tube used in the measure- 
ments just described except that the envelope 


*was of glass of high ultraviolet transmission, in- 


stead of quartz. The average value of r calculated 
from measurements with this tube was double 
that determined with the first tube for the same 
vapor density and electrode voltages. The meas- 
urements were inconsistent and not reproducible. 
With a peak alternating voltage of 4.5 volts on 
G—S and H—P, and temperatures between 200° 


* P. H. Garrett, reference 7, p. 784. 
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and 214°C, the calculated lifetime varied between 
6.4 and 3.8X10~* sec. On the other hand, when 
the peak voltage was decreased to 1.5 volts, the 
calculated lifetime was 2.3 X 10~* sec. Absorption 
measurements showed that the radiation meas- 
ured was still substantially 3261A. These spuri- 
ous results were later found to be due to a 
defective coating on the hot cathode. The greater 
part of the electron emission was from the lower 
tip of the cathode, which was not opposite the 
window of the oven. The bulk of the radiation 
which reached the photoelectric tube was excited 
in this region and could reach the photoelectric 
tube only by being absorbed by the vapor in the 
region opposite the window and reradiated. As a 
result the calculated lifetime was approximately 
twice the true lifetime. The variability in the 
results and their dependence upon voltage was 
due to the variation of the distribution of the 
emission over the surface of the cathode with 
time and applied voltage. When the window was 
lowered to expose the active part of the cathode, 
the measurements gave results in close agreement 
with those found with the quartz tube. 


2288A 


A series of measurements was made to study 
the behavior of the principal series resonance line 
2288A, with the same quartz excitation and 
photoelectric tubes as were used with 3261A. 

The concentration of cadmium vapor in the 
excitation tube was varied between 1.6X10"° 
(113°C) and 2.7 X10" (184°C), atoms cm~*. The 
more precise and significant measurements were 
made at temperatures of the reservoir V below 
140°C and with crest values of the voltage G—S 
less than 10.5 volts. Under these conditions the 
only radiation emitted by the excitation tube, to 
which the photoelectric tube responded, was 
2288A. To test the radiation for the presence of 


3261A, a Corex D filter which absorbs all wave- 


lengths shorter than 2500A was interposed be- 
tween the tubes. When the temperature of the 
reservoir was below 170°C and the crest value of 
G—S was less than 12 volts, no radiation was 
transmitted which could be detected by the 
photoelectric tube; ie., no radiation between 
2500A and 4000A was present in a measureable 
amount. To test for radiation shorter than 2500A, 
other than 2288A, absorption measurements were 


made with the cadmium vapor absorption cell] 
between the excitation and photoelectric tubes, 
The radiation tested was from the excitation tube 
when the reservoir was at 145°C. The tubes were 
separated so that the distance between excitation 
tube and absorption cell was 72 mm and that 
between absorption cell and photoelectric tube 
was 90 mm. With high pressures in the absorption 
cell, absorption was practically complete, indi- 
cating that the radiation affecting the photo- 
electric tube was essentially 2288A. 

With lower pressures in the absorption cell, the 
absorption coefficient of cadmium vapor for this 
radiation was computed to be 1.75X10-"N, 
which agrees, within the precision of the observa- 
tions, with the value, 1.64X10-"N, given by 
Zemansky.® In these calculations allowance was 
made for that part of the radiation scattered 
from the absorption cell which reached the 
photoelectric tube. 

Throughout this series of measurements the 
voltages W—S and K—P were held at 4.5 volts. 
The static characteristics were determined in the 
manner described above. That of the excitation 
tube (radiation as a function of G—S) was very 
nearly a straight line with an intercept on the 
voltage axis at 6.0 volts. The accepted excitation 
potential for 2288A is 5.4 volts. The characteristic 
of the photoelectric tube was of the usual form 
with saturation for H—P less than 1 volt. The 
current reversal occurred at —0.7 volt. 

The R—f curves were made with fixed bias 
voltages; G—S equal to 4.5 and H—P to —0.7 
volt. The alternating voltage impressed upon the 
directly connected external terminals of G and H 
was maintained for all frequencies at 5.0 volts, 
peak value. The reference frequency was 14,000 
cycles. It was found, however, that the radiation 
from the excitation tube decreased with increase 
in frequency as determined by applying alter- 
nating voltages to G—S only, H—P being held 
constant. To correct for this, the ratio of the 
radiation for each frequency to that for the refer- 
ence frequency was determined. The observed 
value of R was then corrected to correspond to 
equal radiation excited at all frequencies, by 
dividing by this ratio. The cause of this variation 
was presumed to be voltage drop in the leads 
within the excitation tube at high frequencies, 
resulting in a voltage difference between the 
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electrodes less than that between their external 
terminals. A decrease in the electrode voltage 
G—S of 3 percent would have accounted for the 
maximum observed variation of the radiation 
emitted, which was about 10 percent. A corre- 
sponding decrease in the voltage on the electrodes 
of the photoelectric tube would have had negli- 
gible effect on the measurements. 

The results of the measurements are shown in 
Fig. 3, which gives, as a function of the concen- 
tration of cadmium atoms, the values of the 
“apparent lifetime’’ of the radiation emitted by 
the excitation tube. In an earlier paper the au- 
thors have shown’ that in the case of radiation 
excited by electron impact, provided (a) that 
absorption is negligible, (b) that the excitation of 
radiation in the excitation tube is proportional to 
the instantaneous value of the alternating voltage 
during the positive half-cycle, and zero during 
the negative half-cycle, and (c) that the current 
in the photoelectric tube passes abruptly from a 
positive saturation value to a negative saturation 
value upon reversal of the alternating voltage, 
the relation between R and f is given by 


(1) 


Here s is the ratio of the negative to the positive 
photoelectric saturation current, and 7 is the 
mean lifetime of the radiation from the excited 
atom. Except for the effect of absorption and 
reradiation at the higher concentrations, these 
conditions were well satisfied in the measurements 
with 2288A as seen from the static characteristics 
and the bias voltages selected for the electrodes. 
When no complications due to absorption are 
involved, r can easily be found from the curves. 
From (1) we find: 


for R=R.=1—}(1+5); (2) 

for 

(3) 
r= (4) 


where f, is the frequency for which R is the mean 
of unity and the final asymptotic level R., ap- 
proached by the curve at high frequencies. In 
Fig. 3, the values of (2xf;)—' were found from the 
experimental curves, and from the value of s 
calculated from the photoelectric characteristic. 


Hence 
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The use of this value instead of the experimentally 
determined value was necessary because of the 
difficulty of making measurements at sufficiently 
high frequencies and because the curves departed 
markedly from the curve as given by the equa- 
tion. The experimental curves dropped more 
sharply than the theoretical curve, and R fell toa 
value considerably lower than R,, as calculated 
from the static characteristic. Although it was 
not possible to verify this by measurements at 
sufficiently high frequency, there were indications 
that this low value of R was a minimum and that 
the curves would have approached the theoretical 
value of R, had higher frequencies been used. 
The extrapolation of the curve in Fig. 3 to zero 
concentration, to eliminate the effect of absorp- 
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Fic. 3. Curve showing persistence of 2288A as a function 
of the concentration of cadmium atoms. The value of 
(2xf\)~ at zero concentration equals 


tion and reradiation, gave 11.0 10~-* sec. for the 
limiting value of (2xf,)—'. This would give a life- 
time of 2288A about five and one-half times as 
great as the value 1.99 X 10~* sec. found by earlier 
investigators.** Because of this large difference 
in the results and because of the unusual form of 
the experimental curves, it seemed probable that 
the process under investigation was more com- 
plex than the simple excitation of 2288A postu- 
lated in the computations. This was confirmed by 
a study of the form of the curve in Fig. 3. For 
brevity, let (2xf,;)-' be represented by T. The 
increase in T with the vapor concentration JN is 
a complex function involving the geometry of the 
tube and cannot be readily computed. The earlier 
work of the authors and that of Garrett? with the 
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2537A line of mercury showed that the relation 
between J and N is given approximately by the 


equation 
T=T)+ANbr, (5) 


where A is the atomic absorption coefficient and } 
the “effective thickness” of the vapor traversed 
by the radiation. In the case of the mercury 
radiation, the constant T> is identical with 7, the 
lifetime of the resonance radiation. From (5) we 
have 

r=(1/Ab)(dT/dN). (6) 


The value of 5 can in general be only roughly 
estimated from the geometry of the tube. How- 
ever, the tube used in this work was nearly 
identical in design and dimensions with that used 
by Garrett, whose results are well described by 
(5) if b is taken as equal to 0.5 cm. Using this 
value of band Zemansky’s value of A, 1.6410-", 
and dT/dN, 1.73X10-*°, determined from the 
curve, we find from (6) the value r=2.110-° 
sec., which, considering the approximations used, 
agrees well with the value 1.99X10~® sec. de- 
termined by other methods. 

The possibility that the peculiar shape of the 
curve was due to a phase shift, causing the 
voltage on the excitation tube grid to lag behind 
the voltage applied to the photoelectric tube 
grid, was investigated. Such a phase shift would 
have resulted in a curve having a final level below 
the R,, calculated from the static characteristics. 
A comparison of the curves at high and low con- 
centrations of cadmium vapor showed equal 
differences between the minimum values of R and 
R. although the frequencies used in the two 
cases differed by a factor of ten. Since the phase 
shift should have decreased with the frequency, 
it was concluded that it was not sufficient to 
cause significant error. 

The most consistent explanation of these re- 
sults is that under the conditions of the experi- 
ment, the excitation of the atoms by electron 
impact was mainly to a level (possibly several 
levels) a few tenths of a volt above the 5'P; level, 
and that excitation of the 5'P; level was by an 
infra-red transition to it from the higher level. 
The persistence of 2288A observed by the method 
used in these measurements would depend, there- 
fore, upon both the lifetime of the 5'P; state and 
on 7; the lifetime of the infra-red transition. Since 


the photoelectric surface was sensitive only to the 
2288A radiation, under the conditions used in 
these measurements the R—f curves for such a 
cascade process are described by the equation 


4x’ tT 4x f?r? 


From this equation it can be shown that with an 
error of less than two percent, To, the value of T 
at zero concentration, determined as described 
above is equal to r+7;. Using the earlier de- 
termined value of 7, 1.99 10~-® sec., we find for 
the lifetime of the infra-red transition, or average 
lifetime of the several transitions if there were 
more than one, 9.010-* sec. The curve repre- 
senting Eq. (7) has the peculiarities observed in 
the experimental curves. It drops more sharply 


than the curve described by (1), falls to a mini-' 


mum well below R,, and at high frequencies ap- 
proaches it as a limiting value. Substituting in it 
the values of r and 71, found above, we find, for 
zero concentration a minimum value of R at 
f=5fi, 0.06 below R, which agrees approximately 
with the observed values. This result is similar to 
that found for the corresponding line in the mer- 
cury spectrum 1849A by the authors.* It was 
shown that when excitation was produced by 
electron impact, practically all of the 1849A radi- 
ation emitted by the tube was due to cascading 
from a level about 0.4 volt above the 6'P, level, 
this infra-red transition having a lifetime of 
2.1X10-* sec. This result agreed with that of 
Eldridge’® who had found by other methods that 
for excitation by electron impact, the bulk of 
1849A was excited indirectly through transitions 
from higher levels. 


Subordinate Series Lines 


Studies made to determine the mean lifetimes 
of several of the subordinate series lines of 
cadmium gave results which were not conclusive, 
but are briefly described because of certain inter- 
esting facts brought out by the measurements. 
These measurements were made with the potas- 


* H.W. Webb and H. A. Messenger, reference 7, p. 326. 
J. A. Eldridge, Phys. Rev. 23, 685 (1924). . 
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sium photoelectric tube previously used and 
described by Randall.’ Filters were employed to 
jsolate the lines being studied. To obtain suffi- 
cient energy in these lines it was found necessary 
to use voltages on the grid of the excitation tube 
which resulted in considerable ionization. There 
js consequently some uncertainty as to how 
closely the excitation followed in time the voltage 
impressed upon the grid. Numerous checks indi- 
cated, however, that this was not a serious source 
of error in the frequency range employed. A more 
serious one was the very considerable absorption 
of the radiation by atoms in the 5*P states which 
in most of the measurements tended to obscure 
the effect of the persistence in the higher energy 
states. Figure 4 shows a typical R—f curve taken 


1 1 
2 3 
FREQUENCY x 


Fic. 4. R—f curve for the group of subordinate series 
lines 5086-4800-4678A. Values of fy occur at 45,000 and 
1.6X10® cycles, for the two sections of the curve, re- 
spectively. 


- under these conditions. Initially the curve falls 
off sharply to a value of R well above the final 
level, R.. The frequencies at which this drop 
occurs corresponded to a persistence one to two 
times as great as that which would be due to the 
lifetime of 3261A. Since the filters between the 
tubes did not transmit 3261A but only the sub- 
ordinate series lines under measurement, it was 
concluded that this initial drop was due to ab- 
sorption of the subordinate series lines by atoms 
in the 5*P states and their subsequent re-emis- 
sion. The decrease of R at very much higher 
frequencies will be discussed below. 
Measurements of the mean lifetime of the 6°S; 
level, from which emission of the lines 5086, 4800, 
and 4678A takes place, were made with a special 
excitation tube. This had an envelope of Pyrex 


glass and a cathode similar to that used in the 
other tubes. It was constructed with two coaxial 
electrodes of nickel mesh, G; and Ge, with radii 
4.5 mm and 8.0 mm, respectively, and a coaxial 
wall electrode of sheet nickel W with radius 18 
mm. The radiation was observed parallel to the 
axis through a window at the end of the tube, and 
baffles in the tube were so placed that only radia- 
tion from the region between Gz and W reached 
the photoelectric tube. Corning filters 368 and 
429 cut off all radiation of wave-length less than 
4100A and isolated the subordinate series lines 
5086, 4800, and 4678A as a group. The bias 
voltages used on the electrodes were: W—S=45 
volts; G.—S=15 volts. The alternating voltage, 
of 15 volts peak value, was applied between G, 
and S. The R—f curve shown in Fig. 4 is typical 
of a group taken with a vapor pressure of 
8.4X10-* mm (216°C). 

The initial drop of R was from 1.00 to 0.69, and 
the frequency for which R had the half-value 
(0.85) for this drop was 45,000 cycles, correspond- 
ing to a lifetime of 3.510~*, which is approxi- 
mately double the lifetime of 3261A. The value of 
R remained nearly constant until the frequency 
500,000 cycles was reached, after which it 
dropped gradually to the final level 0.21. The f; 
for this portion of the curve was taken as the 
frequency, 1.6X10° cycles, for which R was the 
mean of the values 0.69 and 0.21. Assuming that 
this second drop in the curve was determined by 
the lifetime of the 6*S, state and that no absorp- 
tion and re-emission of these subordinate series 
lines was present, we find the lifetime of this state 
is approximately 10~’ sec. (2zf,)—. If absorption- 
re-emission processes were involved, this would 
be simply an upper limit which these experiments 
set for this lifetime. Some measurements made on 
these same lines with the tube of glass of high 
transmission in the ultraviolet, described earlier, 
gave values of f, in the neighborhood of 8X 10° 
cycles, which can be interpreted as the result of 
sufficient absorption and re-emission of the radia- 
tion to increase the average persistence of the 
6°S; state to double that of the true mean life. 

With a set of filters isolating the group of 
diffuse subordinate series lines 5*P — 5*D, similar 
measurements were carried out. The curves are 
in most respects similar to those just described, 
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and show that the value for the lifetime is not 
greater than 10~’ sec. 

A detailed explanation of these curves is diffi- 
cult. The initial drop seems best explained as 
dependent on the subordinate series radiation 
which was excited more or less continuously in 
the region screened from the photoelectric tube 
by the baffles and which was absorbed and 
reradiated by the excited atoms in the 5°P states 
found in the unscreened part of the tube. The 
corresponding persistence would then be the 
same as the persistence of the 5*P states in that 


which occurred at high frequencies, would be 
explained in the case of the sharp series as due to 


part of the tube. The second drop in the curve, 


TABIN 


the direct excitation to the 6°S level (or 5*D leve} 
for the diffuse series) of atoms in the unscreened 
region during the positive half-cycle of the 
alternating voltage and their radiation to the 
photoelectric tube either directly or with not 
more than one or two absorption-re-emission 
steps. For this radiation the persistence would be 
of the order of the mean lifetime of the level in 
question. Because of the uncertainties as to the 
exact processes involved in the measurements, 
the above results can probably be relied upon 
only as to the order of magnitude of the lifetimes, 

The authors take this opportunity to thank 
Dr. Herbert N. Otis for his assistance in this 
investigation. | 
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various producing materials. 


INTRODUCTION 


N 1938 Schein and*Wilson! published a paper 
concerning the production of penetrating 
cosmic-ray secondaries in the atmosphere above 
20,000 feet. They found considerable production 
at 25,000 feet and indicated that the secondary 
particles observed were produced by photons. 
Since then, many attempts have been made to 
identify the nature of the producing radiation. 
Shonka? found the production to be as high as 6.1 
percent in a producing layer of 20 cm of lead at 
an altitude of 14,200 feet. Only 1.5 percent was 
attributed to penetrating radiation arising from 
photons. The other 4.6 percent was ascribed to 
neutral particles of greater penetration (neu- 


1M. Schein and Wilson, Phys. Rev. 54, 304 (1938). 
2 F. Shonka, Phys. Rev. 55, 24 (1939). | 
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Production of Single Mesotrons by Non-Ionizing Radiation at 
Altitudes of 10,600 ft. and 14,200 ft. 
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The production of single mesotrons in various materials was studied by means of an apparatus 
consisting of 46 Geiger-Mueller counters arranged in 10 different three-, four-, and fivefold 
coincidence sets. The producing layers consisted of different thicknesses of paraffin, iron, and 
lead. The results indicate that photons are the most probable agent for the production of the 
single mesotrons. The cross section for the mesotron production has been calculated for the 


trons). Schein, Jesse, and Wollan* undertook the 
investigation of ‘the production of penetrating 
particles by non-ionizing radiation as a function 
of altitude. Their results indicated that an ap- 
preciable production takes place above 23,000 
feet with increasing altitude. Because the in- 
tensity of the produced mesotrons followed the 
intensity of the soft component, photons were 
considered as the most probable producing agent. 
Rossi and Regener‘ studied the production of 
mesotrons at 14,200 feet. They found a definite 
production which was attributed to either 
neutrettos or high energy neutrons. They arrived 
at this conclusion because it appeared that the 
producing radiation had a large probability of 


3 Schein, Jesse, and Wollan, Phys. Rev. 56, 613 (1939). 
* B. Rossi and V. Regener, Phys. Rev. 58, 837 (1940). 
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Fic. 1. Counter arrangement. 


traversing several centimeters of lead without 
generating secondary ionizing particles. Regener,*® 
in 1943, reported that the production of non- 
shower-producing particles by photons at alti- 
tudes up to 11,500 feet was as high as 10 percent 
of the total penetrating component. All these 
different investigations definitely showed that 
there exists a production of mesotrons by non- 
ionizing radiation. However, no definite evidence 
was available to indicate the nature of the 
mesotron-producing radiation. 


’ The following experiment was carried out at 


Mt. Evans, Colorado (14,200 feet), in order to 
determine whether photons are ‘able to produce 
single mesotrons. This seemed possible by in- 
vestigating the absorption of the mesotron pro- 
ducing radiation as a function of thickness in 
materials of different atomic numbers. 


APPARATUS 


The counter arrangement used in the following 
experiment is schematically represented in Fig. 1. 
The forty-six counters were arranged into ten 
different three-, four-, and fivefold coincidence 


*V. Regener, Phys. Rev. 64, 252 (1943). 
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circuits. The seventeen counters A and also the 
eleven counters B were each connected in parallel. 
Counter batteries A and B completely covered 
and extended in both directions beyond the solid 
angle subtended by ceunters CDEF. The coinci- 
dence circuits were of a conventional type similar 
to those used by Schein, Jesse, and Wollan in 
their balloon flights. The resolving power of the 


circuits was tested and found to be approximately .- 


10-5 sec. Hence, for the three-, four-, and five- 
fold coincidences chance coincidences were negli- 
gible. Each coincidence circuit operated an indi- 
vidual neon light which was photographed by a 
Simplex movie camera whenever a pulse occurred 
simultaneously in the so-called master coinci- 
dence set. Counters CDEF constituted’ this 
master set. A coincidence pulse in the master set 
corresponded to the passage of a penetrating 
particle® through the lead layers X. Thus, by 
observing the neon lights which were photo- 
graphed during each exposure on the same film, 
one could determine which of the various possible 
coincidence sets were tripped simultaneously 
with the master set. 

All counters were 1 in. in diameter, but of two 
different lengths. The counter batteries A and B 
consisted of counters 13.5 in. long. All other 
counters were 6 in. long. The tubes were filled 
with a mixture of 90 percent argon and 10 percent 
alcohol. Each counter tube had a voltage plateau 
of at least 200 volts, and all were matched so that 
they had very nearly the same starting potential. 
They all were then connected to the same high 
voltage supply and were operated at 1240 volts. 
The counters were mounted in a wooden stand 
which contained various lead layers X and Y’ and 
provided space for the insertion of a producing 
layer. In Fig. 1 the producing layer is designated 
by 2. This layer could be varied up to a thickness 
of 11.5 cm. = consisted of paraffin, iron, and lead. 
The three lead layers X were used to discriminate 
between penetrating particles and electrons, and 
each had a thickness of either 0.64 cm or 1.59 cm 
of lead. Since a thickness of 0.6 cm is very nearly 
equal to a radiation unit in lead, an electron 
passing through all three lead absorbers should 
have an extremely high probability of emerging 


*Since most of these particles are later shown to be 
mesotrons, they will be referred to as such in the rest of 


* this paper. 
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as a cascade shower below each of the lead layers 
X. To determine whether the particle which 
tripped the master set was an electron or a 
mesotron, several groups of counters were ar- 
ranged in seven separate coincidence sets so as to 
indicate the presence of showers below the vari- 
ous lead layers. Counters GCHF formed a four- 
fold and MNF a threefold coincidence set. Pulses 
in these sets showed the presence of a shower 
directly below the producing layer >. Counters 
IDJF, KELF, and CFST constituted three four- 
fold (coincidence sets), and counters OPF and 
QRF two threefold coincidence sets. Coincidence 
pulses from any of these counter arrangements 
indicated that showers were present below the 
lead layers X. It is to be noted, however, that a 
small percentage of the threefolds occurring 
below only one single layer X could be due to the 
presence of knock-on electrons kicked out of one 
of the lead layers by a mesotron.’ The lead layer 
>’, 6.98 cm thick, was used to absorb low energy 
mesotrons. 

ACDEF formed a fivefold coincidence set. If 
this set was tripped, the presence of an ionizing 
particle traversing the producing layer = and the 
layers X was indicated. A coincidence in CDEF 
which was not accompanied by a coincidence in 
ACDEF generally meant the production of an 
ionizing particle in layer 2 by non-ionizing radia- 
tion. These produced particles were attributed to 
single mesotrons when such an event was unac- 
companied by a pulse in any of the seven shower 
sets. The lower set of counters B was also con- 
nected with the master set CDEF in a fivefold 
coincidence. A coincidence in this set showed 
whether or not a particle tripping the master set 
had enough energy to traverse the 6.98 cm of the 
lead absorber 2’. 

With the six side-counters in bank A removed, 
some cases were found where ionizing particles 
entered from the side and were scattered in the 
producing layer in such a way that they tripped 
the master set. Such events appeared similar to 
the production of an ionizing particle in 2. After 
adding the six “‘police’’ counters in A, however, 
the presence of side particles was detected as a 
coincidence A CDEF and hence did not represent 
a production process by non-ionizing radiation. 


7 This effect proved to be negligibly small for the in- 
frequent number of produced mesotrons. 
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- entire course of the experiment. 


Air showers had no effect on the production of 
single mesotrons because they were eliminated by 
tripping at least one of the shower sets in the 
neighborhood of counters CDEF. The efficiency 
of the banks of counters A and B was very nearly 
100 percent. This high efficiency could be ob. 
tained only by using tubes of relatively small 
diameter, thereby reducing the single counting 
rate in each of the tubes in banks A and B. 


EXPERIMENT 


In order to determine if single mesotrons could 
be produced by non-ionizing rays, over 20,000 
pictures of mesotrons tripping the master set 
were taken at Mt. Evans. From preliminary tests 
conducted at Echo Lake, Colorado (10,600 feet), 
it was determined that the frequency of the above 
process, if it existed at all, was a very small per- 
centage of the frequency of the penetrating 
component. It was necessary to keep all the 
counters in a fixed position during the entire 
course of the experiment. Therefore, in the final 
arrangement, the only variables were the thick- 
ness of the lead layers X and the thickness and 
material of the producing layer 2. The inefficiency 
of the upper bank of counters was tested by re- 
moving the producing layer 2. When this was 
done, all the cases where only the coincidence set 
CDEF was tripped and which were unaccom- 
panied by a simultaneous discharge in set 
ACDEF could then be attributed to inefficiency 
of A. Likewise, the inefficiency of the counter 
bank B was tested by removing lead layer »’, 
These inefficiency tests were repeated several 
times during the course of the experiment and 
were found to average 0.24 percent of the fre- 
quency of the penetrating component.® Since the 
position of the counters was not changed during 
the entire course of the experiment, the actual 
production of single mesotrons could be obtained 
by correcting the measured frequency of produc- 
tion processes for the inefficiency in A and B. 


RESULTS AND INTERPRETATION 


The results on production of single penetrating 
particles on Mt. Evans and Echo Lake are given 
in Table |. A slow particle refers to a mesotron or 
any other penetrating particle of lower energy 


8 The inefficiency did not exceed 0.4 percent during the 
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TABLE I, Data obtained at Mount Evans and at Echo Lake. 


At Mount Evans 


Material Material used for the 


ucing layer 
ieee ‘ a 1 2 3 4* 5 6 
Thickness 0 producing 
1.27 Ae 11.40 11.40 11.40 1.27 


layer = in cm 
1.59 1.59 1.59 0.64 1.59 


layer X in cm 
Number of mesotrons in- 

0107 99.0 95.2 88.6 66.3 
0.61 O41 066 0.74 0.05 


cident on = per hour 
Produced penetrating 
06 0.32 O42 0.19 0.41 


mesotrons per hou 
Produced slow mesotrons 
per hour** 


Lead 
7 8 9 11 


3.81 1140 0.64 1.27 11.40 
1.59 159 064 0.64 0.64 


79.2 66.0 53.3 87.1 684 
0.39 O11 0.09 0.25 0.12 
0.22 0.23 063 1.64 0.37 


——— 


* An additional 57.2 cm of paraffin was placed above counter bank A 


** Inefficiency has been subtracted from the apparent values of the production. 


which trips set CDEF but is unable to traverse 
the absorbing layer >’. The experiment also gave 
interesting information about multiple produc- 
tion in different materials and cascade showers in 
lead. The results of these investigations will be 
reported at a later date. 

The single particles which are produced and 
which traverse the master set CDEF are pene- 
trating particles, providing none of the shower 
sets is simultaneously tripped. Most of these 
penetrating particles must be mesotrons since 
protons would require an average energy of 
4X 10° ev to traverse the lead layers interposed 
between the counters, and there is no evidence 
for the existence of protons of such high energies 
at altitudes of Mt. Evans. Many cloud-chamber 
photographs have been taken by various ob- 
servers, but only protons of lower energy (E < 10° 
ev) were found.’ 

There are different lines of evidence which tend 
to throw light on the nature of the producing 
radiation. (1) Cosmic-ray photons are mainly 
absorbed by pair-production processes with a 
cross section which is proportional to the square 
of the atomic number Z. One should expect that 
if the mesotron-producing radiation consists of 
photons, it would be absorbed according to this 
law. Since the absorption of photons in lead 
(Z=82) is very great, one would expect the pro- 
duction by photons to reach saturation in a lead 
thickness of 2~—3 cm. In paraffin, where Z is very 
much smaller, one would not expect saturation 
until a producing layer of about a meter is 


*In the table a production of particles always refers to 
a production process by non-ionizing radiation. 

D. and S. Neddermeyer, Phys. Rev. 50, 
263 (1936). W Bostick, Phys. Rev. 61, «d (1942). 
W. M. Powell, ‘Phys Rev. 61, 670 (1942). W. E. Hazen, 
Phys. Rev. 65, 67 (1944). ; 


reached. One can see from Table | that in lead, 
and particularly in that part dealing with the 
production of slower mesotrons passing through 
the 0.64-cm lead layers X (columns 9, 10, 11), 
this effect is very well pronounced. The apparent 
decrease in production with 11.4 cm of lead is due 
to the self-absorption of the produced mesotrons 
in the 11.4 cm of lead. In 11.4 cm of paraffin this 
effect does not occur, and one sees that saturation 
is not reached in this thickness. In order to find 
out whether the production is saturated in a still 
larger thickness, an additional paraffin layer of 


‘57.2 cm was placed above the upper counter bank 


A. This paraffin completely covered the solid 
angle of the four counters CDEF of the master 
set. If saturation were reached, most of the pro- 
duced mesotrons would originate in the 57.2-cm 
paraffin block, trip the counter bank A, and 
thereby enter the producing layer of 11.4 cm as 
an ionizing particle. One would, therefore, expect 
a very abrupt drop in the registered production 
of mesotrons by non-ionizing radiation. Re- 
ferring again to Table I, one sees very distinctly, 


Tasve II. Number and en of mesotrons 


produced in lea 


ingle mesotrons 
uced per hr. 


1.64 
0.41 
0.22 
0.37 
0.23 


after comparing the total production of mesotrons 
in column 3 and 4, that this does not occur, 
showing that in paraffin there is no saturation 
effect with a thickness of 57.2 cm. 

(2) The energy range of the mesotrons pro- 
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duced in lead which are stopped in 2’ are com- 
puted from the ionization loss in lead and are 
given in Table II. (The results were calculated 
from the number of produced mesotrons which 
were stopped in the lead layer 2’). In this 
computation it was assumed that all the produced 
mesotrons originated from a layer one cm below 
the top of the producing lead block 2. Table II 
shows that there is a considerable drop in the 
number of produced mesotrons at energies lower 
than 1.6X 10° ev. It should be noted that in the 
photon spectrum at 10,000 and 14,000 feet alti- 
-tudes a similar change occurs around energies of 
about 2X10® ev.'' Assuming a mesotron rest 
mass of 200 electron masses, it takes at least 
10* ev to create a single mesotron. If the proba- 
bility of production does not change appreciably 
with photon energy in the neighborhood of 108 ev, 
there should then be a drop in the energy spec- 
trum of the produced mesotrons at an energy of 
about 10° ev such as observed, provided the 
produced mesotrons are created by photons. 

(3) The increase in production of single meso- 
trons between the altitudes of Echo Lake and 
Mt. Evans can be seen by comparing the data at’ 
Echo Lake (given in columns 14 and 15 of 
Table 1), with the corresponding production at 
Mt. Evans (given in columns 3 and 8). This in- 
crease in the case of paraffin” amounts to a ratio 
of 2.6. The increase in production is roughly of 
the same order of magnitude as the increase in 
the intensity of the soft component (photons) be- 
tween the two altitudes. This ratio is about 2 
according to Greisen." 

Comparing the number of produced single 
mesotrons in lead found in this experiment with 
that obtained under very similar conditions in a 
stratosphere balloon flight carried out in collabo- 
ration with Dr. Schein in 1942, we find at 5 cm 
Hg pressure that the number of mesotrons pro- 
duced in 2 cm of Pb is about 100 times the value 
(Table I, columns 6 and 7) obtained at Mt. 
Evans. In the same balloon experiment the in- 
crease in the photon intensity for energies 


" W. E. Hazen, Phys. Rev. 65, 67 (1944). 
12 A similar comparison in lead cannot be made because 


' the self-absorption in the producing lead layer overshadows 


the major production of slow mesotrons. 

1K. Greisen, Phys. Rev. 61, 212 (1942). 

“The photon intensity was put proportional to the 
measured number of cascade showers produced in 2 cm 
of Pb by non-ionizing rays. 
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higher than 10 ev was measured and found to be 
600 between sea level and altitudes corresponding 
to 5 cm Hg pressure. Since the photon intensity 
increases by a factor of eight between sea level 
and Mt. Evans," one should, therefore, expect an 
increase in the production of mesotrons by non- 
ionizing radiation between Mt. Evans and an 
atmospheric height equivalent to 5 cm Hg to be 
600/8 = 75, which is in as good agreement as can 
be expected under the conditions imposed.'* This 
comparison is very rough and represents only an 
order of magnitude consideration. 

The arguments presented under 2 and 3, while 
not conclusive, definitely favor photons as possi- 
ble producing agents of slower mesotrons. Argu- 
ment 1 substantiates this conclusion strongly 
since photons are the only non-ionizing rays 
which are able to explain the observed differences 
in absorption between paraffin and lead. 

Assuming photons as the producing agent, one 
can calculate the cross sections for the mesotron 
production for the various producing materials, 
From the data obtained in paraffin, it is con- 
sidered probable that the mesotrons are created 
by collisions with the atomic nucleus of the 
producing material. Under these assumptions the 
cross section in smaller thicknesses of paraffin can 
be calculated from the formula, 


where o represents the cross section for 'mesotron 
production by photons; N, the number of 
photons passing through the thickness ¢ without 
being absorbed by mesotron production ; No, the 
number of incident photons; and n, the number 
of nuclei per cc. Changing to logarithms and 
expanding, we get, 


No—N 
No 


1/,No—N\? 
2\ No 3X No 
Since the number of produced mesotrons (No—) 


is very small in comparison with the number of 
photons No striking the producing layer, the 


No 
not =log = 


% The thickness of the producing layer and the counter 
geometry were not identical in the two cases. 
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higher order terms can be neglected and, 
not™(No—N)/No 


and, therefore, 


For paraffin this formula gives the correct 
value for the cross section for thicknesses in which 
the absorption of photons by pair production is 
negligible. This is very nearly the case for thick- 
nesses of 1.27’cm and 3.81 cm. It was estimated 
that if we take the absorption of photons by pair 
production into account, the cross section would 
change by about 5 percent. In lead and iron the 
absorption of photons by pair production is ap- 

iable even in the smallest thicknesses used. 
In this case the formula gives a value which may 
be considered a lower limit for the cross section 
for mesotron production by photons. 

For monatomic substances n= (6.02 10" 
Xdensity)/atomic weight. This gives for lead 
n= 3.2 X10" and for iron n=6.2 X10” nuclei per 
cc. In order to get m for paraffin we use the 
formula C,Hen+42. This leads to a value for n of 
1.0X 10” nuclei per cc. 

From Hazen’s experiments! the number of 
photons with energies higher than 10° ev is about 
10 percent of the total number of mesotrons 
passing through 4.8 cm of lead. The value of No 
was obtained from the measured number, of 
mesotrons incident on the producing layer Y and 
passing through 4.8 cm of lead. 

Using the experimental values for No and N, 


we get for the cross sections per atomic nucleus 
for this process : 


Opar. =5X10-*%% = 1K 107" cm?; 


These values are about equal to the areas of 
the corresponding nuclei calculated by using the 
approximate formula for the radius of a nucleus 
r=roZ',!® where ro is the range of the nuclear 
forces= 1.4 10-". The areas calculated in this 
way are, 


A = 2 X cm?; 
A iron = 1X 107-75 cm?. 


A par. = 2X 10-*5 cm?; 


This result favors the hypothesis that the meso- 
tron production by photons has to be considered 
as a direct interaction between the photon and 
the atomic nucleus. 

In conclusion the author wishes to express his 
appreciation to Dr. Marcel Schein for suggesting 
the problem and for his numerous suggestions 
and consultations throughout this work. He also 
wishes to express his thanks to General L. A. 
Lawson and Captain Innes-Taylor for their 
generous hospitality at Echo Lake, to Mr. 
George E. Cramner and Dr. William Hyslop for 
making possible the stay at Mt. Evans, and to 
Dr. A. H. Compton for his continuous support in 
these investigations. 


‘6 This formula can be considered only as a fair + 
mation for elements of higher atomic numbers, like Pb. 


ct an 
non- 
d an 
to be 
can 
This 
ly an 
yossi- 
\rgu- 
ngly 
rays 
Nces 
, one 
tron 
rials, 
con- 
ated 
: 
s the 
| can 
xi- 
tron 
hout 
, the 
nber 
and 
- N) 
r of 
inter | 


PHYSICAL REVIEW 


Letters to the Editor 


ROMPT publication of brief reports of important dis- 

coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


Gases Evolved by Magnetized Iron 
in Sulfuric Acid 
R. L. Horr, J. J. NAuGHTON, —4 SMOLUCHOWSKI, 
AND H,. H. 


Research Laboratory, General Electric Aenea Schenectady, New York 
May 29, 1944 


EVERAL months ago we were asked to report whether 
oxygen was one of the products resulting from the 
reaction of magnetized iron with dilute sulfuric acid. The 
first experiment was set up as shown in Fig. 1. A piece of 
Swedish iron 4.7 cmX1.2 cm diameter was placed in 5 
percent sulfuric acid within the field of a solenoid whose 
strength at the center measured 286 oersteds. The acid 
was prepared from distilled water boiled vigorously for 
more than one hour and cooled in nitrogen to remove 
appreciable dissolved oxygen. Gas evolved by the iron 
specimen was collected in a glass-stoppered tube and 
sealed to a microgas-analysis equipment designed for the 
analysis of gases in metals. Oxygen was determined by 
change in pressure of the dried gas, after reacting with 


Fic. 1. Experimental arrange- 
ment to collect gases from mag- 
— iron in dilute sulfuric 
a 
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collection of gases from pole 
pieces P. 


Fic. 2. Arrangement for the (a) 


SCALE 


hydrogen over hot platinum and condensing out water in a 
solid CO: trap. The complete analysis was as follows. 


Oxygen 0.0+0.2 vol. percent 
Hydrogen 94.7 vol. percent 
Hydrocarbon as methane 4.5 vol. percent 
Residue as nitrogen 0.5 vol. percent 


The method was sensitive to better than 0.2 vol. percent 
oxygen, but within this limit, none was found. 

Another experiment was performed in the apparatus 
shown in Fig. 2. The field strength used was now much 
higher, namely, 13,000 oersteds (measured in air) between 
two 0.64-cm (} in.) diameter pole pieces P of Swedish iron 
separated by a 0.16-cm (7g in.) gap. Water likewise was 
boiled and cooled in nitrogen, and at no time was it in 
contact with air. Gases were collected in tube B in which 
there was always some water above the stopcocks in order 
to prevent gas leakage. Analyses were performed by two 
methods: One*was the same as that described above and 
gave the same negative result; the second method using 
150-ml samples was that of oxygen absorption in pyrogallol. 
The result of the latter analysis was 0.02 to 0.05 vol. per- 
cent oxygen. The difference between gas collected with and 
without magnetic field was within the limit of error and 


random in sign. 


On the Quantization of Probability 


Henry F. DuNLAP 
University of New Mexico, Albuquerque, New Mexico 
March 2, 1944 
HE recently proposed quantization of probability! 
leads to an apparent contradiction when applied toa 
particular problem. 
Consider an enclosure of volume V, containing a grain 
of sand of volume G. The enclosure is assumed field free 
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except for the gravitational field of the grain of sand. 
If v is any subvolume of V large enough to contain G, then 
the a priori chance of finding the grain of sand at any 
instant in v is 

P=0/V.. 

Now suppose that V is made so large that P is less than 
one chance quantum. According to the postulate of 
quantized probability, the grain of sand cannot exist in v. 
But v was any subvolume of V, hence the grain of sand 
cannot exist anywhere in V. An attempt to remove the 
contradiction by assuming that geometrical space is 
quantized in such a way that the smallest allowable volume 
js greater than v fails, since v was initially taken to be larger 
than G, the volume of the grain of sand. 


1 Alfred N. Goldsmith, Phys. Rev. 64, 377 (1943). 


Comments on “On the Quantization of 
Probability” by Henry F. Dunlap 


ALFRED N. GOLDSMITH 
580 Fifth Avenue, New York, New York 
March 12, 1944 


HILE the writer would prefer to leave the analysis 

of such thoughtful comments relative to the conse- 

quences of the probability-quantization theory as are given 

in the foregoing commentary to active specialists in the 

field of probability theory, he desires nevertheless to offer 
the following brief suggestions. 

In the commentary a conclusion is reached that, under 
the stated conditions and in accord with probability 
quantization, the hypothetical “grain of sand cannot 
exist” any longer. At first sight, this seems a common- 
sense conclusion derivable from the premises. 

But if probability is indeed quantized, ‘‘common sense” 
in its everyday meaning will not serve in such an analysis, 
any more than ‘common sense” particularly and reliably 
assists us in the realm of quantum physics. It is possible to 
examine the problem from a different viewpoint and to 
draw a different conclusion. Let us for example first sub- 
stitute for ‘‘grain of sand” the probably more appropriate 
term ‘‘smallest elementary particle.” Such a particle might 
be of electronic dimensions or even less. Let us also substi- 
tute for ‘‘cannot exist” the phrase “‘cannot be located"’ or 
“cannot be discerned.’’ This substitute language seems 
more in accord with the necessary treatment of the 
problem. We can then paraphrase the above commentary 
by stating that quantized probability appears to require 
that a sufficiently small particle enclosed in an adequately 
large volume becomes unlocatable by any physical means. 
This conclusion is not found troublesome or self-contra- 
dictory by this writer. 

It may seem a radical or strange conclusion to state 
that when an extremely small particle is enclosed in a 
sufficiently large space it is, in a physical sense, altogether 
lost according to our chance-quantum theory. This would 
be tantamount to saying that it has become literally un- 
findable by any means available to the would-be observer. 
Yet this conclusion seems no more radical than others 
freely accepted by modern physicists. Our “lost particle” 
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would be no more unfindable than is an amount of energy 
less than one quantum or an amount of matter less than 
the smallest elementary particle, both according to present- 
day theory. Our “lost particle’ has no presence in our 
physical world in any meaningful sense, and no experiment, 
no matter how intellectually ingenious, will serve to 
locate it. The physical quest for such a particle thus 
becomes as meaningless as the search for an absolute frame 
of reference for motion. Such at least might be the in- 
evitable consequence of probability quantization as the 
writer sees it, and he finds nothing intellectually repugnant 
in such a conclusion. 

It should, however, be stressed, as it was in the original 
letter describing probability quantization, that no assertion 
is made as to the correctness and validity of the theory of 
quantized chance. It was there urged that this theory 
continue to be thoroughly explored, keeping in mind, 
however, that the final conclusion will doubtless be reached 
by experimentation. 

The above commentary is appreciated by the writer 
and is believed by him to be a stimulus to further study of 
the subject. It is hoped that such study will deal with 
methods for the experimental proof or disproof of proba- 
bility quantization as well as with attempted applications 
of such a theory to certain small-scale or infrequent elec- 
trical phenomena. In these last-mentioned realms the new 
theory might display special utility or significance. 


Magnetic Ions 


V. D. Hopper 
Department of Physics, University of Melbourne, Melbourne, Australia 
March 27, 1944 


HRENHAFT?! has described experiments which he 

claims establish the existence of “magnetic ions.” 
His method involves a measurement of the velocity of small 
particles (radii approximately 10~* cm) of various elements 
(Ni, Fe, Mn, Cr, Sb) which are acted on by the earth’s 
gravitational field and a uniform vertical magnetic field 
produced by an electromagnet. The particles are illumi- 
nated at right angles to the direction of the field and are 
observed in a direction normal to both the field and the 
direction of illumination. 

An attempt has been made to verify his conclusions by 
slightly modifying the apparatus used by T. H. Laby and 
the author* in the determination of the electronic charge, 
a horizontal magnetic field being substituted for the 
electric field used in that determination... The magnetic 
field was produced by a small permanent Alnico magnet 
which has a strong field strength of about 1000 oersteds 
at the midpoint between the poles. Fine particles of nickel 
(radii of the order of 10-* cm) were drawn by an air current 
between the poles of the magnet and then allowed to settle 
under the action of the vertical gravitational and the 
horizontal magnetic fields, and their path was observed 
by means of a microscope. Thousands of particles have 
been observed, many being in the field of view at the one 
time, but no particle has been detected moving in such a 
manner as to verify the “magnetic ion” hypothesis. In 
the center of the magnetic field all particles fell vertically, 
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and on either side of the central line particles either re- 
mained undeflected or deflected towards the nearer pole. 
No particles crossed the center of the field. This result 
agrees with existing theories of magnetic induction for the 
non-uniform but symmetrical field between the poles of 
the magnet. The experiment has been repeated with iron 
and rouge particles with the same results. Increase in 
illumination also made no difference except that vertical 
convection currents were produced when excessive heat 
entered the chamber. 

A possible explanation of Ehrenhaft’s results is that the 
particles he observed were electrically charged and were 
acted on by some stray electric field such as that produced 
between the coils of the electromagnet. A typical particle 
observed by Ehrenhaft has a radius of approximately 
10-* cm and maximum velocity due to the electromagnet 
of 10-* cm/sec. Assuming that this velocity is due to an 
electric field and that the particle possesses a charge of ne 
where n is the number of electrons each of charge 
e(+4.8X10-" e.s.u.) and applying the corrected form of 
Stokes’s law, one finds that the field required to produce 
this velocity is of the order of 0.3/n e.s.u. or 90/n volt/cm. 

Thus for sufficiently highly charged particles a velocity 
of the above magnitude would be observed for very small 
electric fields (of say <1 volt/cm) and these small fields 
may be produced by the potential differences between the 
two halves of the electromagnet or electrostatic leakage 
effects from the generator. 


1 Felix Ehrenhaft, Phys. Rev. 37, 659 (1940); Ann. d. Physik lle, 
Sec. 13, 151 (1939-40) ; ce 96, 228 (1942); J. Frank. Inst. 253, 225 


(1942). 
?V. D. Hopper and T. H. Laby, Proc. Roy. Soc. A178, 242 (1941). 


Test for Change of Pole Strength of 
Permanent Magnet 


J. E. GotpMAN 
Westinghouse Research Laboratories, East Pitisburgh, Pennsylvania 
May 22, 1944 


T the Pittsburgh meeting of the American Physical 
Society, F. Ehrenhaft' reported an approximately 10 
percent loss of strength of a permanent magnet in an experi- 
ment in which acidulated water is placed between the 
poles of the magnet in contact with the pole faces and 
allowed to remain so for a period of time. The energy loss 
of the magnet is claimed by Ehrenhaft to have gone into 
the “magnetolysis” of the water. In view of the unusual 
character of these results, it was decided to repeat Ehren- 
haft’s experiment as described in the Bulletin of the A meri- 
can Physical Society. 

In our experiment, a stabilized Alnico permanent magnet 
was used, the strength of which was measured on five 
successive days before the commencement of the experi- 
ment in the same laboratory where the entire experiment 
was conducted. Measurements were made by means of a 
ballistic galvanometer and search coil, the ballistic galva- 
nometer having been calibrated by a standard of mutual 
inductance prior to each measurement. The accuracy of 
the apparatus is better than 1 percent. A special attach- 
ment was constructed to assure that in each measurement 
the search coil would start and come to rest at precisely 
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the identical position. Pole pieces for the Alnico magnet 
were made of Armco iron, the end of each terminating in 
a truncated cone the face of which was 6 mm in diameter, 
The geometrical arrangement is shown in Fig. 1. Measure. 
ments of the field before and after the experiment were 
carried out both with and without the pole pieces. The 
height of each pole piece was carefully measured under a 
microscope by comparison with a standard height gauge, 
The accuracy of this measurement was +0.004 mm. The 
purpose of this last phase of the investigation was to follow 
any possible change in gap length that might accompany 
the action of the acid on the soft iron pole faces. 

Two experiments were carried out to test directly Ehren- 
haft’s results. In one, 4 percent H2SO, was placed in the 
glass cell, and the iron pole faces were exposed to the acid, 
As expected, there was a violent display of bubbles due to 
the liberation of hydrogen by the iron. At the end of 18 
hours, the field of the magnet was measured and found to 
be identical with the value obtained at the beginning of 
the experiment within the measure of experimental error. 

In the second experiment, the conditions were the same 
except that a 12 percent sulfuric acid solution was used in 
place of the weaker one used previously. This concentration 
was estimated to be of approximately maximum electro- 
lytic activity. The entire pole structure with all but the 
pole face covered with paraffin was immersed in the solu- 
tion and allowed to remain so for a similar period of time. 
No measurable decrease in field was observed as a result of 
this experiment. It is estimated that from the commence- 
ment of the investigation until its completion, the magnet 
poles were exposed to the “magnetolytic” action of the 
acid for approximately 60 hours during all of which time 
the magnet was not removed from its original surroundings, 
i.e., from the same laboratory where the measurements 
were made. The total magnetic flux at the pole face of the 
original Alnico magnet was 11,700+50 maxwells. This 
value was the same at the end of the experiment as at the 
beginning. It may be added parenthetically that when the 
experiment was carried out with the pole faces covered with 
paraffin and the acid previously boiled to drive off dissolved 
gases, no bubbles at all were observed when the acid was 
placed in the cell whereas they are observable if the acid 
is not previously boiled. 

1F, Ehrenhaft, Phys. Rev. 65, 349 (1944). 
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